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A COMPUTER PROGRAM FOR THE DETERMINATION OF THE 
ACOUSTIC PRESSURE SIGNATURE OF HELICOPTER ROTORS 
DUE TO BLADE THICKNESS 

G, H. Mall* and F. Farassat** 

Langley Research Center 

SUMMARY 

This report presents a computer program for the determination of the thickness 
noise of helicopter rotors. The results are obtained in the form of an acoustic pressure 
time history. The parameters of the program are the rotor geometry and the helicopter 
motion descriptors. The formulation employed is valid in the near and far fields. The 
blade planform must be rectangular, but the helicopter motion is arbitrary. The observer 
position is fixed with respect to the ground with a maximum elevation of 45° above or 
below the rotor plane. With these restrictions, the program can also be used for the cal- 
culation of thickness noise of propellers. 

INTRODUCTION 


Thickness noise of helicopter rotors or propellers is the noise generated by the 
normal velocity distribution on the surface of these rotating bodies. In the past, thick- 
ness noise of helicopter rotors and propellers has not been studied thoroughly even after 
the introduction of some approximations. This noise was investigated by Deming (ref. 1) 
and Arnoldi (ref. 2); the latter based his analysis on the work of Billing (ref. 3). Deming's 
and Arnoldi '3 analyses contain some of the following restrictions which should be removed 
in applications. These restrictions are of three types: 

(1) The rotor or the propeller is required to be stationary or to move in axial direc- 
tion at uniform speed. 

(2) The observer is located in the far field. 

(3) The acoustic source distribution must be compact. 
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For example, Deming's analysis is developed with restrictions 1 and 2. The expression 
derived in references 4 and 5 removes restrictions 1, 2, and 3. This paper documents a 
computer program developed to evaluate this expression. 

The program, as presented, is for blades with rectangular planform. A variety of 
uniform airfoil sections can be defined by the user in the program input data. The 
observer is fixed with respect to the ground; this condition is preferred when experi- 
mental and theoretical results are compared. Numerical examples showing favorable 
comparison with some experimental results are presented in reference 5. 


THE FORMULATION 


Let the surface of the rotor system be specified by f(Y,r) = 0, where Y is a 
Cartesian frame fixed to the undisturbed medium and t is the time. If v^ is the local 
normal velocity of the blade surface and p is the acoustic pressure, then the governing 
equation for the determination of the thickness noise is (refs. 4 and 5) 


c^ 9 t2 






( 1 ) 


where c is the speed of sound in the undisturbed medium, 6 (f) is the Dirac delta func- 
tion, and Pq is the density of the undisturbed medium. The solution of equation (1) is 
presented in reference 4 and in slightly more general form in reference 5. 

In the present work, the 77 '- and ^-frames in reference 5 are replaced by the rf- and 
YH-frames, respectively. Consider one single blade and a nonrotating coordinate system: 
the YH-frame is fixed to the center of rotation (fig. 1), and a rotating system is the 
?f-frame, fixed to the blade (fig. 2). Let T(t 7 j,t 72 ) be the thickness function of the blade 

in this frame. Then the solution of equation (1) is given by equation (32) of reference 5 
as follows: 


p(X,t) 


8 c ^2 r 

277 9t Jr(D.p.) 


TlVi + T2V2 



dr dT 


( 2 ) 


where 

p(X,t) acoustic pressure 

Pq density of undisturbed medium 
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x,t 


observer position and time 


Tj,T 2 source times when the sphere g 

the blade, respectively (X and 

Y source location Y-frame 

(VHi,VH 2) vehicle velocity components along Y^ and Y2 axes of Y-frame 
n angular velocity of rotor 

(fi,f2,r3) unit vector along radiation direction in if- frame 

curve of intersection of sphere g = 0 and mean surface of the blade 
disk plane 

= dT/dVv T2 = aT/9r?2 
= -VHi + r?2fi 
= -VHg - 

For several blades, the contribution of each to the acoustic pressure is added linearly. 

THE PROGRAM 
General Discussion 

For this program, the blades are assumed to have a rectangular planform with a 
uniform airfoil section along the span. The thickness distribution function T(rjj,?72) is 
denoted as F(t7i) here since there is no dependence on 772- The Rotor Thickness Noise 
program (RTN) evaluates the integral in equation (2) numerically. Evaluation of this double 
integral involves one integration over the curve of intersection of the sphere g = 0 and 
the rotor blade; the other integration is over source time. The first of these integrals is 
evaluated using a trapezoidal integration. The second integral is performed with Simpson's 
rule and a trapezoidal end-point correction where required. The Rotor Thickness Noise 
program postprocessor (RTNPP) completes the calculation of the acoustic pressure by 
differentiating the RTN output numerically with respect to observer time. The purpose of 
this section is to present detailed descriptions of the RTN and RTNPP programs, to 


r 

D.P. 

Tl 

Vl 

^2 


= T- t + |x-Y|yc = 0 enters and leaves 
t fixed), |x - y| = r 
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describe the manner in which input data must be prepared, and to describe the program 
output. Flow charts and FORTRAN listings of the programs are included as appendixes. 


The solution technique employed by the RTN program is presented in the section of 
this report entitled "Method of Solution.'/ Within this section, several key internal pro- 
gram variables are identified parenthetically. The RTN program reads a number of 
parameters describing the rotor system structure, the amount of detail to be used in per- 
forming the integrations, and various program options. A detailed description of the input 
required by the program is given in the section entitled "Preparation of Input Data." Out- 
put from the RTN program includes an echo of the input data and a table of ordered pairs 


[q5(x,t),t 


where 


is the observer time and T’(X,t) is defined by equation (3). 


27t Jtj ‘^r(D.P.) 


ViF'ivi) 


1-v^K F'2(t,i)(i - f ^2 


1/2 


dr dT 


( 3 ) 


The section entitled "Postprocessor Description" discusses the numerical differentiation 
methodology used by the RTNPP program in computing the thickness noise from 
equation (3). 


Method of Solution 

Based upon the translational velocity of the rotor system (VH) and the specified 
initial value of source time (TAU), the location of the center of the rotor system relative 
to a fixed coordinate system is determined. This fixed coordinate system (Y) is described 
in figure 1. The origin of this system is the location of the center of the rotor system at 
time TAU = 0. The Yj^Y2-plane contains the rotor disk and the Yj-axis is parallel to the 
initial orientation of blade number one. The direction from the blade tip to the hub is the 
positive direction. A new reference frame (YH) is constructed with its origin coinciding 
with the center of the rotor system at time TAU. This nonrotating coordinate system 
remains fixed to the helicopter. The position of the observer at time TAU in the YH-frame 
(X) is computed from the specified location in the Y-frame (XO). The blade-fixed coordi- 
nate system is shown in figure 2. 

At time TAU, a sphere is constructed with its center at the position of the observer 
with a radius (RM) equal to the maximum distance between the observer and any blade tip. 
The observer time (OTIME) is determined by the relation OTIME = TAU - RM/SNDSPD, 
where SNDSPD is the speed of sound. The radius of this sphere is subsequently permitted 
to collapse at the speed of sound. The points of intersection of the projection of this 
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sphere on the rotor system plane and the circle defined by the outer radii of the blades 
are determined as shown in figure 3. At the position of the observer, the more clockwise 
of these points (YSTART) defines the initial sweep angle (PHIO). The other point of inter- 
section (YEND) defines the final sweep angle. 

The integral over the intersection of the sphere and the rotor blade is performed by 
sweeping the arc of the intersection of the projection of the sphere and the rotor system 
plane. The arc is swept in a counterclockwise sense. To execute the sweep of such an 
arc efficiently, the rotor system is divided into four regions as illustrated in figure 4. 
These regions are defined as follows: 

IREGION = 0 The point lies outside the circle defined by the outer radius of the 
blades. 

= 1 The point lies between the circles defined by the outer and inner 
radii of the blades, but is not on a blade. 

= 3 The point lies inside the circle defined by the inner (hub) radius 
of the blades. 

= 6 The point lies on a blade. 

Three sweep increments are used by the program. In region 1, a rather coarse 
angular increment (DELPHIl) is used to space across the region. When a blade is 
encountered, the intersection of the arc with the edge of the blade is computed. A 
smaller angular increment (DELPHI2) is used to sweep across the blade. An even 
smaller angular increment is allowed for points near the leading edge for blunt leading 
edges. The values of the integrand and differential increment are calculated for each 
point on a blade. The line integral (trapezoidal) is computed dynamically by cumulative 
summing. The individual contributions of each blade are also summed. The hub region 
is bypassed by using the computed value of the intersection points of the arc with the 
inner circle. The line integral is completed when the sweep across the arc reaches the 
final sweep angle. 

When a sweep has been completed, the value of the sweep time (SWPTAU) is 
increased by a small increment (DELTAU). The center position of the rotor system ' 
relative to the observer is recomputed. A new sphere is constructed and its line integral 
is calculated. This process is continued until the sphere no longer intersects the rotor 
system. As these line integrals are calculated, they are dynamically accumulated into a 
Simpson's rule sum. The integrands are multiplied by two or four if they are even or odd 
terms, respectively, and are then added to the sum. If the last term is even, a trapezoidal 
correction is made. 

After the last line integral has been calculated, the value of the double integral and 
the corresponding observer time are output. The value of the source time is increased 
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by some multiple of DELTAU and the entire process is repeated until one of three possible 
termination criteria is satisfied. These criteria are: 

(1) TAU exceeds a-specified m^imum value. 

(2) The radius of the contracting sphere becomes less than the outer radius of the 

blades. ^ 

(3) The angle of the observer relative to the center of the rotor system and the 
rotor system plane exceeds 45°. 

The hierarchy of program RTN and its 10 subprograms is illustrated in figure 5. Flow 
charts and FORTRAN listings of every RTN subprogram are presented in appendix A. 

Preparation of Input Data 

All input data to RTN are associated with the single NAMELIST name ROTOR. 

These variables are summarized in table I according to dimension, type, and units and 
are described in the following paragraphs. 

CH Blade chord. 

COEFFS The coefficients for describing the airfoil shape. The airfoil shape 

used in the subprogram FOFETAl is given in nondimensional form by 

F(E)/CH = T(Ci\fE + C2E + C3E2 + C4E3 + CsE^) 

where E = t/j is the distance along the chord from the leading 

edge, CH is the blade chord, and T is the thickness ratio. The function 
F(E) is the airfoil thickness function which is the airfoil ordinate at t]^. 
This equation is suitable for many airfoil shapes, in particular for the 
NACA four digit airfoil series. The coefficients for this series (from 
ref. 6) are as follows: 

Cl = 1.4845 

C2 = -0.6300 

C3 = -1.7580 

C 4 = 1.4215 

C5 = -0.5075 

DELTET The angular displacement of the rotor system between consecutive posi- 

tions of the contracting sphere. This variable is used with the rotor 
system angular velocity to compute DELTAU, the time required for the 
rotor to rotate DELTET degrees. DELTAU is then used as the differ- 
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DTAUM 


ETAMAX 


JFLG 


KFLG 


MORDATA 


N 


ential increment in the integration over time. The source time TAU is 
increased by some multiple of DELTAU between computations of 
cp(X,t). A useful formula to select DELTET is 

DELTET = (CHMOMEGA) cos » 

2.5(345 - VH cos \p) 


where ify is the angle that the line from observer to the center of rota- 
tion makes with the vertical and VH is the helicopter speed. 

A factor to determine the time increment between data points. The 
time increment derived from the variable DELTET is multiplied by 
DTAUM to determine the time interval between data points. Conse- 
quently, the rotor system rotates DTAUM+DELTET degrees between 
data points. A useful rule for selecting this factor is 


DTAUM = 


4 to 6 

(NBLADES)(DELTET) 


The maximum value for which increased resolution is desired 
near the leading edge for airfoils with blunt leading edges. Generally, 
a value between 10 to 15 percent of the chord is sufficient. 

Rotor system translational motion indicator: 

JFLG = 0 Hovering helicopter 

= 1 Helicopter in forward flight 

Airfoil section type indicator: 

KFLG = 0 Sharp leading edge 
= 1 Blunt leading edge 

Additional data indicator: 

MORDATA = .T. Process another data set 
= .F. Terminate 

A factor to control the spacing along the arc for points not on a blade 
(i.e., to skip from one blade to the next). The variable N should not be 
made too small since a large sweep increment near a blade tip might 
skip over the blade completely. Approximately N points are used to 
traverse a distance along the arc equal to the blade chord. A value 
between 3 and 8 is suggested. 
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NB LADES 
NS 


OMEGA 

PERDM 


R 

RO 

SNDSPD 

TAUINT 


THKRAT 

XO 


The number of blades in the rotor system. The program allows a max- 
imum of eight blades. 

A factor to control spacing along the arc for points on a blade. The 
variable NS is used to compute the sweep angular increment, DELPHI2. 
The DELPHI2 determines the line integral differential increment. 
Approximately NS points are used to traverse a distance aloi^ the arc 
equal to the blade chord. A value between 15 and 25 is suggested. 

Angular velocity of the rotor system. 

The number of periods of the blade rotation to be processed. Because 
of the lengthy computation time, this number should not be given a large 
value. This variable is used to establish the maximum time termina- 
tion criterion. 

Radius of the blades. 

Hub radius. 

The speed of sound. 

The initial value of TAU, the source time. Nonzero values of this vari- 
able provide a restart capability as well as a means of generating high- 
resolution calculations in TAU regions of interest. 

Blade thickness ratio. 

The location of the observer relative to (0,0,0), the initial position of the 
center of the rotor system at time TAU = 0 (i.e., the Y-frame). 


Postprocessor Description 

The relationship between the RTN program and the postprocessor (RTNPP) is illus- 
trated in figure 6. Sample RTN output which is input to RTNPP is shown in figure 7. No 
input, other than the information generated on TAPE 12 by RTN, is required by the 
postprocessor. 

The Rotor Thickness Noise program postprocessor in use at Langley Research 
Center (LRC) performs three functions. The RTN output ordered pairs (<T ,t) are used to 
find the thickness noise by differentiating them numerically with respect to t, the observer 
time. The resulting thickness noise is interpolated to provide a specification at equally 
spaced observer time values. The results are then plotted. 

Because cp is the result of numerical integration, error resulting from finite resolu- 
tion should be eliminated before attempting to differentiate the RTN output. This elimi- 
nation is accomplished with a subroutine which fits a smooth cubic spline curve to the set 


8 



« 


of data points. This routine requires a specification of the allowed variation of each 
functional value. For each data point, the difference between the raw data and a value 
obtained by trigonometric smoothing is used to satisfy this requirement. 

The postprocessor output consists of four plots (illustrated in figs. 8 to 11) for each 
data set. The raw data as output by RTN are displayed in figure 8. The cubic spline 
curve is presented in figure 9. Figure 10 displays the smoothed data over the raw data 
to provide a measure of the quality of the cubic spline fit. Figure 11 shows the acoustic 
pressure time history. 

Flow charts and FORTRAN listings of several RTNPP subprograms are presented 
in appendix B. Several routines required by RTNPP were obtained from the LRC math 
and graphics libraries. Descriptions of these subroutines are included in appendix B. 

CONCLUDING REMARKS 

This report presents a computer program for calculating the thickness noise of 
helicopter rotors with rectangular blades. Favorable comparisons with experimental 
results containing high-speed helicopter blade slap have been found and are reported in 
reference 5. The program is written in Control Data Corporation (CDC) FORTRAN and, 
although several postprocessor subroutines are not fully documented because of propri- 
etary rights, sufficient information is supplied to adapt the program to any facility. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
December 12, 1975 
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APPENDIX A 


ROTOR THICKNESS NOISE PROGRAM 

This appendix provides a brief description of the calculations performed by each 
RTN subprogram. The subprogram flow chart and FORTRAN listing follow the 
description. 


Program RTN 

RTN reads and writes the NAMELIST input data and initializes other variables used 
elsewhere in the program. Following each contraction of the sphere, a check is made to 
see if any termination criteria have been satisfied or if the sphere has completed its pass 
through the rotor system. If the sphere is still intersecting the rotor system, an addi- 
tional contribution is made to the Simpson’s rule integral. If the sphere has passed out 
of the rotor system, the source time is incremented, the computed value of the double 
integral is output, and an additional sphere contraction is initiated. If one of the termina- 
tion criteria has been satisfied, the program terminates unless an additional data set is 
specified. 
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APPENDIX A 


Listing of Program RTN 


c* 

c* 

c» 

c* 

c* 

c* 

c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c» 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c* 
c * 
c* 


PROGRAM RTN(OUTPUT=191B,INPUT=l01B,TAPEl2=OUTPUT,TAPE5=INPUn 


ROTOR THICKNESS NOISE 


PURPOSEI 

PROGRAM RTN OETERMINES THE ACOUSTIC PRESSURE 
GENERATED BY A ROTATING THICK BODY, SUCH AS AN 
IDEALIZED RIGID ROTOR, HITH SOME SPECIFIED MOTION 
RELATIVE TO AN OBSERVER AT SOME INITIAL DISTANCE 
FROM THE BODY. THE PRESENT PROGRAM INCLUDES 
ROUTINES DEVELOPED SOLELY FOR ROTOR SYSTEMS KITH 
RECTANGULAR BLADES AND UNIFORM AIRFOIL SECTION 
ACROSS THE SPAN OF THE BLADES. 


METHOD 

OF 

SOLUTIONI 

THE FORM OF THE SOLUTION UTILIZED BY THIS PROGRAM 
INCLUDES A DOUBLE INTEGRAL WHERE ONE INTEGRATION 
IS OVER THE INTERSECTION OF THE CONTRACTING 
SPHERE C*<TAU-T) *R=0 (CENTERED AT OBSERVER) 

AND THE BODY, AND THE OTHER INTEGRATION IS OVER 
SOURCE TIME, TAU. THE FIRST OF THESE INTEGRALS IS 
EVALUATED USING A TRAPEZOIDAL INTEGRATION. THE 
SECONO INTEGRAL IS PERFOPMEO WITH SIMPSONS 
RULE AND A TRAPEZOIDAL END POINT CORRECTION 
WHERE REQUIRED. 

BEGINNING WITH AN INITIAL VALUE OF TIME 
(TAU), THE TRANSLATIONAL VELOCITY OF THE ROTOR 
SYSTEM IS USED TO LOCATE THE POSITION OF ITS 
CENTER. THE LOCATION OF THE OBSERVER RELATIVE 
TO THIS CENTER POSITION IS DETERMINED AND A 
SPHERE IS CONSTRUCTED (CENTERED AT THE OBSERVER) 
HITH A RADIUS EQUAL TO THE MAXIMUM DISTANCE 
BETWEEN THE OBSERVER AND ANY BLADE TIP. THIS 
TIME REPRESENTS THE LOWER LIMIT ON THE TIME 
INTEGRAL AND THIS SPHERE REPRESENTS THE INITIAL 
POSITION OF THE CONTRACTING SPHERE. THE POINTS 
OF INTERSECTION OF THE PROJECTION OF THIS 
SPHERE ON THE ROTOR SYSTEM PLANE AND THE CIRCLE 
DEFINED BY THE BLADES OUTER RADIUS ARE DETERMINED. 
AT THE OBSERVERS POSITION THE MORE CLOCKWISE OF 
THESE POINTS DEFINES THE INITIAL SHEEP ANGLE. 

THE ARC THROUGH THE ROTOR SYSTEM PLANE COR- 
RESPONDING TO THIS POSITION OF THE SPHERE IS 
SWEPT IN A COUNTERCLOCKWISE SENSE TO COMPUTE 
THE LINE INTEGRAL FOR THIS VALUE OF SOURCE TIME. 

TO DECREASE EXECUTION TIME, TWO SWEEP INCREMENTS 
ARE USED BY THE PROGRAM. FOR POINTS ON THE ARC 
WHICH ARE NOT ON A BLADE A LARGE ANGULAR INCREMENT 
IS USED TO SPACE ACROSS THESE REGIONS. SMALLER 
ANGULAR INCREMENTS ARE AUTOMATICALLY USED BY THE 
PROGRAM WHEN A BLADE IS ENCOUNTERED. FOR BLUNT 
EDGES A smaller INCREMENT MAY BE GENERATED FOR 
POINTS BETWEEN THE LEADING EDGE AND A MAXIMUM VALUE 
OF THE DISTANCE FROM THE LEADING EDGE. FOR EACH 
POINT ON A BLADE THE VALUE OF THE INTEGRAND AND 
DIFFERENTIAL INCREMENT APE COMPUTED. THE LINE 
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c* 


INTEGRAL (TRAPEZOIDALI IS COMPUTEO DYNAMICALLY 

• RTN 

67 

c* 


BY CUMULATIVE SUMMING. WHEN A POINT ON THE ARC 

• RTN 

68 

c* 


IS PEACHED WHICH IS OUTSIDE THE CIRCLE DEFINED 

* RTN 

69 

c* 


BY THE ROTOR SYSTEM, THE VALUE OF THE RESULTANT 

* RTN 

70 

c* 


INTEGRAL IS MULTIPLIED BY THE APPROPRIATE PHYSICAL 

* RTN 

71 

c* 


CONSTANT. 

• RTN 

72 

c* 



* RTN 

73 

c» 


THE VALUE OF SOURCE TIME IS INCREASED SLIGHTLY. 

* RTN 

74 

c* 


THE CENTER POSITION OF THE ROTOR SYSTEM RELATIVE 

* RTN 

75 

c* 


TO THE OBSERVER IS RECOMPUTED. A NEW SPHERE IS 

• RTN 

76 

c* 


CONSTRUCTED AND THIS LINE INTEGRAL IS CALCULATED. 

• RTN 

77 

c* 


THIS PROCESS IS CONTINUED UNTIL THE PROJECTION 

* RTN 

78 

c* 


OF the sphere has COMPLETELY PASSED THROUGH THE 

* RTN 

79 

c* 


ROTOR SYSTEM. THE FINAL VALUE OF SOURCE TIME 

* RTN 

8C 

c* 


REPRESENTS THE UPOER LIMIT ON THE TIME INTEGRAL. 

• RTN 

81 

c* 



• RTN 

82 

c* 


SINCE THE SOURCE TIME DIFFERENTIAL INCREMENT IS 

* RTN 

83 

c* 


CONSTANT, A SIMPSON PULE INTEGRATION IS POSSIBLE. 

• RTN 

84 

c* 


THIS INTEGRAL IS PERFORMED DYNAMICALLY AS A 

• RTN 

85 

c* 


CUMULATIVE SUM. AS LINE INTEGRALS ARE CALCULATED 

» RTN 

86 

c* 


THEY ARE MULTIPLIED BY TWO(FOUR) IF THEY ARE 

• RTN 

87 

c* 


EVEN(ODO) TERMS AND ADDED TO THE SUM. IF ThE LAST 

• RTN 

88 

c* 


TERM IS EVEN, A TRAPEZOIDAL CORRECTION IS MADE. 

* °TN 

89 

c* 



* RTN 

90 

c* 


THE VALUE OF TAU IS INCREASED AND THE ENTIRE 

* RTN 

91 

c* 


PROCESS IS REPEATED UNTIL ONE OF THREE POSSIBLE 

• RTN 

92 

c* 


TERMINATION CRITERIA IS SATISFIED. THESE CRITERIA 

* PTN 

93 

r* 


APE 1 

* RTN 

94 

c* 



• RTN 

95 

c* 


1-TAU EXCEEDS A SPECIFIED MAXIMUM VALUE 

• RTN 

96 

c* 



• RTN- 

97 

c* 


2-THE RADIUS OF THE CONTRACTING SPHERE 

• RTN 

98 

c* 


BECOMES LESS THAN THE RADIUS OF THE BLADES 

* RTN 

99 

c* 


CUTER RADIUS 

• RTN 

100 

c* 



• RTN 

ICl 

c* 


3-THE ANGLE OF THE OBSERVER RELATIVE TO 

* RTN 

1C2 

c* 


THE CENTER OF THE ROTOR SYSTEM AND THE 

* RTN 

103 

c* 


ROTOR SYSTEM PLANE EXCEEDS 45 DEGREES 

* RTN 

1'4 

c* 



• RTN 

1C5 

c* 


ONCE ONE OF THE TERMINATION CRITERIA IS SATISFIED 

• RTN 

ICE 

c* 


THE PROGRAM TERMINATES UNLESS AN ADDITIONAL DATA 

* RTN 

10 7 

c* 


SET IS SPECIFIED, NUMERICAL DIFFERENTIATION OF THE 

* RTN 

1L8 

c* 


OUTPUT, AS WELL AS PLOTTING THE RESULTS AND 

* RTN 

1C 9 

c* 


FORMATTING THE RESULTS FOR SUBSEQUENT ANALYSIS 

* RTN 

lie 

c* 


IS PERFORMED IN A COMPANION POST-PROCESSOR. 

• RTN 

111 

c* 



• RTN 

112 

c* 



• PTN 

113 

c* 

PSRAMETERSI 


• RTN 

114 

c* 



» RTN 

115 

c* 

xo 

THE LOCATION OF THE OBSERVER RELATIVE TO (0,0,0), 

* RTN 

116 

c* 


THE initial POSITION OF THE ROTOR SYSTEM AT TIME 

• RTN 

117 

c* 


TAU = 0. (METERS) 

* PTN 

118 

c* 



• RTN 

119 

c* 

NBLADES 

the NUMBER OF BLADES IN THE ROTOR SYSTEM. THE 

* RTN 

12C 

c* 


PROGRAM ALLOWS A MAXIMUM OF 8 BLADES. 

• RTN 

121 

c* 



* RTN 

122 

c* 

R 

RADIUS OF THE BLADES. (METERS) 

* RTN 

123 

c* 



• RTN 

124 

c* 

CH 

BLADE CHORD. (METERS) 

* RTN 

125 

c* 



* RTN 

126 

c* 

OMEGA 

ANGULAR VELOCITY OF THE ROTO® SYSTEM. (REVOLUTIONS 

* RTN 

127 

c* 


PER MINUTE) 

* RTN 

128 

c* 



• RTN 

129 

c* 

RO 

HUB RADIUS. (METERS) 

* RTN 

130 

c* 



* RTN 

131 

c* 

THKRAT 

blade THICKNESS RATIO. 

• RTN 

132 

c* 



• RTN 

133 
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c» 

QELTET 

ANGULAR DISPLACEMENT OF THE ROTOR SYSTEM BETWEEN 

» RTN 

134 

c* 


CONSECUTIVE POSITIONS OF THE CONTRACTING SPHERE. 

• RTN 

135 

c* 


THIS VAEIABLF IS USED WITH THE ROTOR SYSTEM ANGULAR 

♦ RTN 

136 

c» 


VELOCITY TO COMPUTE THE SOURCE TIME DIFFERENTIAL 

* RTN 

137 

c* 


INCREMENT, DELTAU. (DEGREESI 

» RTN 

138 

c» 



* RTN 

139 

c» 

N 

FACTOR TO CONTROL SPACING ALONG ARC BETWEEN BLADES 

* RTN 

140 

c» 


(I.E.. TO SKIP FROM ONE BLADE TO THE NEXT). 

* RTN 

141 

c» 


APPROXIMATELY N POINTS ARE USED TO TRAVERSE A 

* RTN 

142 

c* 


distance along the ARC EQUAL TO THE BLADE CHORD. 

• RTN 

143 

c» 



♦ RTN 

144 

c» 

NS 

FACTOR TO CONTROL SPACING ALONG ARC ON A BLADE. 

♦ RTN 

145 

c» 


APPROXIMATELY NS POINTS ARE USED TO TRAVERSE A 

* RTN 

146 

c* 


DISTANCE along THE ARC EQUAL TO THE BLADE CHORD. 

♦ RTN 

147 

c» 


THIS variable DETERMINES THE LINE INTEGRAL 

» RTN 

148 

c» 


DIFFERENTIAL INCREMENTS. 

• RTN 

149 

c» 



* RTN 

150 

c* 

KFLG 

AIRFOIL SECTION TYPE INDICATOR. 

♦ RTN 

151 

c* 



* RTN 

152 

c» 


KFLG = Qt SHARP LEADING EDGE 

* RTN 

153 

c» 



♦ RTN 

154 

c* 


= 1« BLUNT LEADING EDGE 

* RTN 

155 

c» 



♦ RTN 

156 

c* 

ET&MAX 

MAXIMUM ETA VALUE FOR WHICH INCREASED RESOLUTION 

• RTN 

157 

c» 


IS DESIRED NEAR THE LEADING EDGE FDR BLUNT LEADING 

• RTN 

158 

c* 


EDGES. 

• RTN 

159 

c» 



♦ RTN 

160 

c» 

6LNTN 

FACTOR TO GENERATE FINE SPACING NEAR THE LEADING 

♦ RTN 

161 

c» 


EDGE FOR BLUNT LEADING EDGES. FOR ETA VALUES 

• RTN 

162 

c* 


BETWEEN ZERO AND ETAMAX. THE SPACING FACTOR NS IS 

• RTN 

163 

c» 


MULTIPLIED BY 9LNTN. 

• RTN 

164 

c» 



♦ RTN 

165 

c» 

JFLG 

ROTOR SYSTEM TRANSLATIONAL MOTION INDICATOR. 

» RTN 

166 

c» 



♦ RTN 

167 

c* 


JFLG = Gl HOVERING HELICOPTER. 

• RTN 

168 

c« 



• RTN 

169 

c» 


= 11 HELICOPTER IN FORWARD FLIGHT. 

• RTN 

170 

c» 



♦ RTN 

171 

c* 

COFFFS 

COEFFICIENTS IN THE FUNCTION DESCRIBING THE AIRFOIL 

• RTN 

172 

c* 


SECTION. A MINUS SIGN MUST APPEAR EXPLICITLY FOR 

• RTN 

173 

r.* 


NEGATIVE TERMS. 

• RTN 

174 

c» 



♦ RTN 

175 

c» 

TAUINT 

the' INITIAL VALUE OF SOURCE TIME (TAUI. (SECONDS) 

• RTN 

176 

c» 



♦ RTN 

177 

c* 

PEROM 

THE NUMBER OF PERIODS OF THE ROTOR SYSTEM TO BE 

• RTN • 

178 

c» 


PROCESSED. THIS VARIABLE IS USED TO ESTABLISH 

• RTN 

179 

c* 


THE MAXIMUM TIME TERMINATION CRITERION. 

* RTN 

ISO 

c» 



• RTN 

181 

c» 

DTAUM 

time INCREMENT BETWEEN DATA POINTS. THE SOURCE 

• RTN 

182 

c* 


time differential INCREMENT DERIVED. FROM OELTET 

» RTN 

183 

c» 


IS MULTIPLIED BY DTAUM TO DETERMINE THIS INCREMENT. 

♦ RTN 

184 

c» 



♦ RTN 

185 

c» 

SNHSPD 

SPEED OF SOUND. (METERS/SEC) 

♦ RTN 

186 

e* 



♦ RTN 

187 

c» 

^OROATA 

A logical variable to SPECIFY IF AN ADDITIONAL 

» RTN 

188 

c» 


DATA SET IS TO BE PROCESSED. 

* RTN 

189 

c» 



» RTN 

190 

c» 



♦ RTN 

191 

c» 

REOUIPEO 


* RTN 

192 

c» 

ROUTINES* 


• RTN 

193 

c» 


SUBROUTINE nEWYH 

* RTN 

194 

c* 


SUBROUTINE CALCVH (SEE ADDITIONAL DETAILS) 

♦ RTN 

195 

c* 


SUBROUTINE ORSPOS 

» RTN 

196 

e* 


SUBROUTINE RMOTIME 

* RTN 

197 

c» 


subroutine CALCYST 

♦ RTN 

198 

c* 


SUBROUTINE CALCSP 

♦ RTN 

199 

c» 


FUNCTION FOFETAl (SEE ADDITIONAL DETAILS) 

• RTN 

200 
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FUNCTIOK PFUNCN 


RTN 

201 

c» 


SUaROUIIAE TEST (SEE ADDITIONAL DETAILS) 

♦ 

RTN 

202 

c* 


subroutine LEETA 

* 

RTN 

203 

c» 




RTN 

204 

c» 




RTN 

205 

c* 

AOaiTIONAL 


* 

RTN 

2C6 

c* 

05 T4USI 



RTN 

20 7 

c» 


FOR A NCVING ROTOR SYSTEM (JFLG = 1). A SUBROUTINE 

* 

RTN 

208 

c» 


MUST BE FROVIDEO TO 0ESCRI8E THE MOTION. 

• 

RTN 

209 

c» 




RTN 

210 

c» 


SUBROUTINE CALCVH(VH.TIME) 


RTN 

211 

c* 




RTN 

212 

c* 


SUBROUTINE CALCVH OEFINES THE VECTOR VH(3) 


RTN 

213 

c* 


WHICH MUST contain THE THREE COMPONENTS OF 

« 

RTN 

214 

c* 


THE ROTOR SYSTEM TRANSLATIONAL VELOCITY AT 

9 

RTN 

215 



T = TIME. THESE COMPONENTS ARE TO BE SPECIFIED 

9 

RTN 

216 

c* 


IN THE POTOR SYSTEM REFERENCE FRAME. THE THREE 

9 

RTN 

217 

c» 


AXES IN THIS COORDINATE SYSTEM ARE DEFINED AS 

9 

RTN 

218 

c* 


FCLLOWS , 

9 

RTN 

219 

c* 



9 

RTN 

220 

c» 


THE Y1 AXIS IS PARALLEL TO THE INITIAL 

9 

RTN 

221 

c» 


ORIENTATION OF BLADE NUMBER ONE WITH THE 

9 

RTN 

222 

c» 


DIRECTION FROM THE BLADE TIP TO THE CENTER 

9 

RTN 

223 

c* 


OF THE ROTOR SYSTEM TAKEN AS POSITIVE. 

9 

RTN 

224 

c* 



9 

RTN 

225 

c* 


THE Y2 AXIS IS PERPENDICULAR TO Vl PASSING 

9 

RTN 

226 

c» 


THROUGH THE CENTER OF THE ROTOR SYSTEM AND 

9 

RTN 

227 

c* 


LYING IN THE ROTOR SYSTEM PLANE. THE 

9 

RTN 

228 

c» 


DIRECTION TO The left OF AN OBSERVER FACING 

9 

RTN 

229 

c'» 


rCWAPC POSITIVE ri is taken as positive. 

9 

RTN 

230 

c» 



9 

RTN 

231 

c» 


THE YT AXIS IS PERPENDICULAR TO THE PLANE 

9 

RTN 

232 

c* 


OF THE ROTOR SYSTEM WITH THE POSITIVE 

9 

RTN 

233 

c* 


DIRECTION determined BY THE USUAL CONVENTION 

9 

RTN 

234 

c» 


FOP A right-handed coordinate system. 

9 

RTN 

235 

c* 



9 

RTN 

236 

c» 


THE PRESENT VERSION OF FUNCTION FOFETAl IMPOSES 

9 

RTN 

237 

c* 


A LIMITATION TO UNIFORM AIRFOIL SECTION ACROSS THE 

9 

RTN 

238 

c» 


SPAN OF THE BLADES. FOR NONUNIFORH AIRFOIL SECTION 

9 

RTN 

239 

c* 


THIS function MUST BE MODIFIED. 

9 

RTN 

240 

c» 



9 

RTN 

241 

c» 


THE PRESENT VERSION OF SUBROUTINE TEST IMPOSES A 

9 

RTN 

242 

c* 


LIMITATION TO ROTOR SYSTEMS WITH RECTANGULAR 

9 

RTN 

243 

c* 


BLADES. F0° OTHER BLAOE SHAPES THIS SUBROUTINE 

9 

RTN 

244 

c* 


MUST BE MODIFIED. 

9 

RTN 

245 

c* 



9 

RTN 

246 

c* 



9 

RTN 

247 

c» 

?^ou9ce> 


9 

RTN 

248 

c* 


RTN HAS PROGRAMMED AT LANGLEY RESEARCH CENTER BY 

9 

RTN 

249 

c* 


COMPUTER SCIENCES .CORPORATION BASED UPON A DESIGN 

9 

RTN 

250 

c* 


BY F. FflFASSAT. 

9 

RTN 

251 

c* 



9 

RTN 

252 

c» 



9 

RTN 

253 

c» 

language J 


9 

RTN 

254 

c* 


FORTRAN 

9 

RTN 

255 

c» 



9 

RTN 

2So 

c* 



9 

RTN 

257 

c» 

gate 


9 

RTN 

258 

c* 

RELEASEPI 


9 

RTN 

259 

c* 


may 1. 1974 

9 

RTN 

26D 

c» 



9 

RTN 

261 

c* 



* 

RTN 

262 

c» 

LATEST 


9 

RTN 

263 

c* 

REVISIOMI 


9 

RTN 

264 

c* 


CECEMBEF IE, 1975 

9 

RTN 

265 

c.» 



9 

RTN 

266 

c» 



9 

RTN 

267 

c»* 



99 

RTN 

268 
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APPENDIX A 

LOGICat MORDATA RTN 


DIMENSION 

BETA (8) , 

COEFFS(5) , 

ETA(2,8) , 


RTN 

1 

X(3) . 

XO (3) 

t 

Y(3) , 


RTN 

2 

YENO (2) , 

YH<3) 

« 

YSTART (2) , 


RTN 

3 

VH(3) 





RTN 

COMMON 

ALFA, 

BETA, 

BLADES, 

BLNTN, 

CBFTA, 

RTN 

1 

CH, 

COEFFS, 

OELTAU, 

ETA, 

ETA2, 

RTN 

2 

ETAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

RTN 

3 

IFLAG, 

IREGION, 

J2, 

JFLG, 

KFLG, 

RTN 

4 

NBL AOES, 

OMEGAR, 

OTIME, 

PI, 

R, 

RTN 

5 

PM, 

PMIN, 

S, 

SBETA, 

SNDSPO, 

RTN 

6 

SRCDST, 

SHPCST, 

SHPTAU, 

TAU, 

THKRAT, 

RTN 

7 

TIMTHRU, 

VH, 

X, 

XO, 

YENO, 

RTN 

8 

YH, 

YSTART, 

OELETAl, 

DELF 


RTN 

NAMELIST/ROTOR/XO,N0LAOES, R.CH, OMEGA.RO. THK 

RAT,OELTET 

,N, NS, KFLG, 

RTN 


1 ETAMAX.BLNTN, JFLG.COEFFS.TAUINT.PERDM, DTAUM,SN0SP0,M0R0ATA RTN 

RTN 

READ NAMELIST DATA RTN 

RTN 

10 REAO(5, ROTOR) RTN 

WPITE(12,R0T0R) RTN 

RTN 

INITALIZE variables RTN 

RTN 

PI = 3.1415^265358979 RTN 

OHEGAR = 2.*0MEGA*PI/6C . RTN 

TESTIM = 0. RTN 

YH(1) = 0. RTN 

YW(2) = 0. RTN 

YH<3) = 0. RTN 

Yd) = 0. RTN 

Y(2) = C. RTN 

Y<3) = 0. RTN 

VH<1) = 0. RTN 

VH(2) = 0. RTN 

VH(3) = 0. RTN 

ALFA = ASIN(C.5*CH/R) RTN 

PMIN = R RTN 

RSO=R*R RTN 

POSa=RO*PO RTN 

TAU = TAUINT RTN 

IF(KFLG.EQ.O) BLNTN = 1. RTN 

NPTS=0 RTN 

PKEPl = 0. RTN 

PKEP2 = C. RTN 

SCRIDl = 0. RTN 

SrRIPTP=0.0 RTN 

BLADES = FLOAT(NBLAOES) RTN 

RTN 

COMPUTE SOURCE TIME DIFFERENTIAL INCREMENT AND MAXIMUM TIME RTN 

TERMINATION CRITERION RTN 

RTN 

OELTAU = DELTET/ (6.»0MEGA) RTN 

PERIOD = 6C ./ (0MEGA*9LA0ES) RTN 

TAUMAX = TAUINT ♦ PEROM*PFPIOO RTN 

RTN 

S. N. 2C BEGINS A NEW DATA POINT RTN 

RTN 

20 IF (JFLG.NE. 0) CALL NEHYH(TAU) RTN 

RTN 

DETERMINE RELATIVE LOCATIONS OF POTOR SYSTEM AND OBSERVER RTN 

PTN 

CALL OBSPOS RTN 

PTN 

CONSTRUCT INITIAL SPHERE RTN 

RTN 

CALL RMOTIME RTN 

OTIM = OTIME - TESTIM RTN 

C RTN 
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26 9 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 
28C 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 
295. 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 
3C7 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 
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r 


IF OBSERVER TIME IS PECREaSING SKIP TO NEXT DftT/S POINT 

RTN 

338 

c 



RTN 

339 



IF(DTIM.LE.0.) go to 40 

RTN 

340 



TESTIM = OTIME 

RTN 

341 



OPSftNGL = SSIN(X(3)/SQRT(X(1)*»2 ♦ X(?»**2 * X(3)»*2)) 

RTN 

342 

c 



RTN 

343 

c 


IF ANY TERMINATION CRITERION IS SATISFIED SKIP TO THE 

RTN 

344 

c 


CHECK FOR AOPITIONAL DATA 

RTN 

345 

c 



RTN 

346 



IF (08SANGL. GT. (PI/4. n GO TO 60 

RTN 

347 



IF(RM.LT.RMINI GO TO 6J 

RTN 

348 



IF (TAU.GT.TAUMAX) GO TO 6C 

RTN 

349 

c 



RTN 

35(? 

c 


initialize SMEEP VARIABLES 

RTN 

351 

c 



RTN 

352 



SWPTAU=TAU 

RTN 

353 



TIMTMRU=C.C 

RTN 

354 



IEVEN = -1 

RTN 

355 

c 



RTN 

356 

c 


S. N. 3C begins a new POSITION OF THE SPHERE 

RTN 

357 

c 



RTN 

358 


3C 

TIMTHPU = TIMTHRU ♦ 1. 

RTN 

359 



IEVEN = -IEVEN 

RTN 

360 



FN = N 

RTN 

361 



FNS = NS 

RTN 

36 2 



CALL CALCYST<PSQ,PHIO» 

RTN 

363 



PHI = PHIC 

RTN 

364 

c 



PTN 

365 

c 


ADD TO CUMULATIVE SIMPSON SUM 

RTN 

366 

c 



RTN 

367 



PKEP2 = PKEPl 

RTN 

368 



PKEPl = SCRIPl 

RTN 

369 



PBELT = 2.*SCPIP1 

RTN 

37C 



IF (IEVEN. GT.O ) POELT = 2.*P0ELT 

RTN 

371 



IF(TIMThRU.LT. 2.5» POELT = SCRIPl 

RTN ■ 

372 



SCPIPTP = SCRIPTP ♦ POELT 

PTN 

37 3 

c 



RTN 

374 

c 


IF IFLAG=C, HE ARE IN THE ROTOR SYSTEM 

RTN 

375 

c 



RTN 

376 



IF(IFLAG.EQ.C) go to 50 

RTN 

377 

c 



RTN 

378 

c 


IF IFLAG=2, A termination CRITERION HAS DETECTED IN CALCYST 

RTN 

379 

c 



RTN 

380 



IF(IFLAG.EQ.2» GO TO 60 

RTN 

381 

c 



RTN 

382 

c 


OTHERWISE, HE HAVE COMPLETED THIS DATA POINT. CORRECT THE 

RTN 

383 

c 


SIMPSON integration ANO OUTPUT THE RESULTS 

RTN 

384 




RTN 

385 



POELT = -PKEPl 

RTN 

386 



IFdEVEN.GT.OI POELT = (PKEP2 - 5.*PKEPl»/2. 

RTN . 

387 



SCRIPTP = SCRIPTP ♦ POELT 

RTN 

388 



SCRIPTP = SCRIPTP/3. 

RTN 

389 



NPTS=NPTS+1 

RTN 

390 



P = SCRIPTP*415./(2.*PI) 

RTN 

391 



HRITE(12,70» NPTS,TAU, OTIME, P 

RTN 

392 



SCRIPl = 0. 

RTN 

393 



PKEPl = 0. 

RTN 

394 



PKEP2 =0. 

RTN 

395 



SCRIPTP=0.0 

RTN 

396 

c 



RTN 

397 

c 


INCREASE THE TIME ANO BEGIN A NEW DATA POINT 

RTN 

398 

c 



RTN 

399 


4C 

TAU = TAU ♦ DTAUM*OELTAU 

RTN 

400 



GO TO 20 

RTN 

401 

c 



RTN 

402 

c 


HE ARE IN THE ROTOR SYSTEM. SHEEP THE ARC 

RTN 

403 



■ 
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c 

50 CALL CALCSP<RSO,POSQ,PHI,Y,SCRIPH 
GO TO 30 

CHECK FOR ADDITIONAL DATA 

60 NPTS = 0 

MRITE(12.70I NPTS 
IF(NORDATA) GO TO 10 

TERMINATION. 

70 FORMAT(I5,3E25.15l 
END 


RTN 

4C4 

RTN 

405 

RTN 

406 

RTN 

40 7 

RTN 

4C8 

RTN 

409 

RTN 

410 

RTN 

411 

RTN 

412 

RTN 

413 

RTN 

414 

RTN 

415 

RTN 

416 

RTN 

417 


Subroutine NEWYH 

This routine calculates YH, the three coordinates of the center of the rotor system 
relative to its initial location. This routine is called only for moving (translational) 
rotor systems. 


Flow Chart of Subroutine NEWYH 
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Listing of Subroutine NEWYH 


SUBROUTINE NEHVH(TIME) 


C** 

c» 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c*» 


SUBROUTINE NEHYH IS CALLED BY PROGRAM RTN, AND SUBROUTINE 
CALCYST. THIS ROUTINE RETURNS YH<3». THE THREE 
COORDINATES OF THE CENTER OF THE ROTOR SYSTEM RELATIVE 
TO ITS initial LOCATION (O.OtOI. SUBROUTINE NEHYH IS 
CALLED ONLY FOR MOVING (TRANSLATIONAL* ROTOR SYSTEMS TO 
DETERMINE THE POSITION AT THE TIME SPECIFIED BY THE 
VARIABLE TIME (M/SECl . THE ROTOR SYSTEM VELOCITY (POSSIBLY 
TIME-DEPENDENT* IS OBTAINED FROM SUBROUTINE CALCVH. 


NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 

NEHYH 


a 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


DIMENSION 

BETA (8* , 

COEFFS (5* , 

ETA(2,8* , 


NEHYH 

15 

1 

X(3* , 

X0(3) 




NEHYH 

16 

2 

YE NO (2* , 

YH(3t 

• ' 

YSTART42* , 


NEHYH 

IT 

3 

VH(3I 





NEHYH 

18 

COMMON 

ALFA, 

BETA, 

BLADES, 

BLNTN, 

CBETA, 

NEHYH 

19 

1 

CH, 

COEFFS, 

OELTAU, 

ETA, 

ETA2, 

NEHYH 

2C 

2 

ETAHAX, 

FAC, 

FN, 

FNS, 

FETAl, 

NEHYH 

21 

3 

IFLAG, 

IREGTON, 

J2, 

JFLG, 

KFLG, 

NEHYH 

22 

4 

NBLAOES, 

OMEGAR, 

OTIME, 

PI, 

R, 

NEWYH 

23 

5 

RM, 

RMIN, 

S, 

SBETA, 

SNDSPD, 

NEHYH 

24 

6 

SRCOST, 

SWPDST, 

SHPTAU, 

TAU, 

THKRAT, 

NEHYH 

25 

7 

TIMTHRU, 

VH, 

X, 

XO, 

YENO, 

NEHYH 

26 

8 

YH, 

YSTART, 

OELETAl, 

OELF 


NEHYH 

27 

CALL CALCVH(TIME> 





NEHYH 

28 

DO 10 1=1,3 






NEHYH 

29 

YH(I* = VH(I* 

’TIME 





NEHYH 

3C 

10 CONTINUE 






NEHYH 

31 

RETURN 






NEHYH 

32 

END 






NEHYH 

33 


Subroutine CALCVH 

CALCVH returns the three components of the translational velocity of the rotor sys- 
tem. In general, this subroutine is a user-supplie(i subroutine which must be provided for 
any case other than for that of a hovering helicopter. The routine included with the cur- 
rent version of RTN provides for a constant velocity of 61.73 m/sec (120 knots) directed 
along the Y^-axis. 


Flow Chart of Subroutine CALCVH 
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Listing of Subroutine CALCVH 


SUBROUTINE CALCVH(TIME) 

(^********if*******>t***************9*9‘f*9****************ff****¥*********9* 

C* » 

C* SUBROUTINE CALCVH IS CALLED BY SUBROUTINE NEWYH. CALCVH » 

C* RETURNS VH(3), THE THREE COMPONENTS OF THE VELOCITY * 

C* (TRANSLATIONAL) OF THE ROTOR SYSTEM. THE VARIABLE TIME * 

C* MAY BE USED TO SPECIFY A TIME-OEPENOEN T VELOCITY. * 


CALCVH 

CALCVH 

CALCVH 

CALCVH 

CALCVH 

CALCVH 

CALCVH 


2 

3 

4 

5 

6 

7 

8 


C* * CALCVH 9 

(^¥¥ ¥¥¥¥¥¥*¥¥¥¥»¥**¥¥¥¥¥»*¥¥¥¥ ¥*¥*••¥¥¥*»***¥¥¥¥¥¥¥*¥¥*¥¥¥¥•»»*»¥¥**¥*¥*¥ CALCVH 10 


DIMENSION 

BETA(8), 

COEFFS (5) , 

ETA (2, 8) . 


CALCVH 

11 

1 

X (3) , 

XO ( 3) 

f 



CALCVH 

12 

2 

YENO (2 ) , 

YH(3) 

t 

YSTART (2) , 


CALCVH 

13 

3 

VH (3) 





CALCVH 

14 

COMMON 

ALFA, 

BETA, 

blades. 

BLNTN, 

CBET'A, 

CALCVH 

15 

1 

CH, 

COEFFS, 

DELTAU, 

ETA, 

ETA2, 

CALCVH 

16 

2 

ETAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

CALCVH 

17 

3 

I FLAG, 

I REGION, 

J2, 

JFLG, 

KFLG, 

CALCVH 

18 

4 

NBL ADES, 

OMEGAR, 

OTIME, 

PI, 

R, 

CALCVH 

19 

5 

RM, 

RMIN, 

S, 

SBETA, 

SNOSPD, 

CALCVH 

20 

6 

SRCDST, 

SWPOST, 

SHPTAU, 

TAU, 

THKRAT, 

CALCVH 

21 

7 

TIMTHRU, 

VH, 

X, 

XO, 

YENO, 

CALCVH 

22 

8 

YH, 

YSTAPT, 

OELETAl, 

OELF 


CALCVH 

23 

VH(1) = 61. 

73 





CALCVH 

24 

VH(2) =0. 0 






CALCVH 

25 

VH(3)=0. C 






CALCVH 

26 

RETURN 






CALCVH 

27 

END 






CALCVH 

28 


Subroutine RMOTIME 

This subroutine determines the range from the observer to the most (distant blade 
tip in the rotor system. This range (RM) is the initial radius of the contracting sphere. 
Corresponding to this range and the current value of source time, the value of the 
observer time (OTIME) is also computed. 
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Initialize 
blade counter 
to 1, TOST to 
zero 
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Listing of Subroutine RMOTIME 


c»» 

c* 

c* 

n* 

c* 

c* 

c* 

c» 

c** 


SUBROUTINE RMOTIME 

* ** * * *9 *9 9 * ******** ******¥***»* 9 *** ****** * 

SUBROUTINE RMOTIME IS CALLED BY PROGRAM RTN. TMIS 
DETERMINES THE RANGE TO THE MOST DISTANT TIP IN 
SYSTEM, WHICH IS USED AS THE MAXIMUM SWEEPING 01 
FOR A GIVEN TAU. THE OBSERVERS TIME (CORRESPOND! 
IS ALSO PETUPNEO. 


***•*»•*»•** 

* 

ROUTINE * 

THE POTOR • 
STANCE * 

NG TO TAU) * 


DIMENSION BETA(8), COEFFS 

1 X(3), X0(3), 

2 YEN0(2), YH(3), 

3 VH(3) 


(S>, £TA(2,8), 

YSTART (2» , 




RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 

RMOTIME 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


COMMON ALFA, 

BETA, 

BLAOES, 

BLNTN, 

CBETA, 

RMOTIME 

16 

1 CH, 

COEFFS, 

OELTAU, 

ETA, 

ETA2, 

RMOTIME 

17 

2 ETAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

RMOTIME 

18 

3 IFLAG, 

IREGION, 

J2, 

JFLG, 

KFLG, 

RMOTIME 

19 

4 NBLADES, 

OMEGAR, 

OTIME, 

PI, 

R, 

RMOTIME 

2C 

5 RM, 

rmin. 

S, 

SBETA, 

SNOSPD, 

RMOTIME 

21 

6 SRCOST, 

SWPOST, 

SWPTAU, 

TAU, 

THKRAT, 

RMOTIME 

22 

7 TIMTHRU, 

VH, 

X, 

XC, 

YENO, 

RMOTIME 

23 

8 YH, 

YSTAPT, 

OELETAl, 

OELF 


RMOTIME 

. 24 

DIST(A,B,C,0,E) = (A - C 

l»COS(E))»*2 ♦ (B - 

0*SIN(E) ) 

••2 ♦ C»C 

RMOTIME 

25 

TOST = 0, 





RMOTIME 

26 

XBETA = 0MEGAR»TAU 





RMOTIME 

27 

DO 10 I=1,N3LA0ES 





RMOTIME 

28 

TIPI = XBETA ♦ PI»(1. 

2,*FLOAT<I 

-D/BLAOES) - ALFA 


RMOTIME 

29 

TIP2 = TIPI ♦ 2,»ALFA 





RMOTIME 

30 

n = DIST(X(1) ,X<2) ,X(3) 

,R, TIPI) 




RMOTIME 

31 

IF(DI.GE.TOST) TOST = 01 





RMOTIME 

32 

01 = 0IST(X(1),X<2) ,X(3) 

,R,TIP2) 




RMOTIME 

33 

IF(0I.G£,T0ST) TOST = 01 





RMOTIME 

34 

1C CONTINUE 





RMOTIME 

35 

RM = SORT(TOST) 





RMOTIME 

36 

OTIME = TAU * PM/SNOSPD 





RMOTIME 

37 

RETURN 





RMOTIME 

38 

ENO 





RMOTIME 

39 


Function FOFETAl 

This subprogram computes the blade thickness for specified distance along the 
blade chord. The function included with the current version of RTN is limited to uniform 
airfoil section across the span. 


Flow Chart of Function FOFETAL 
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Listing of Function FOFETAl 



FUNCTION FOFETfiKXETS) 





FOFETAl 

2 

c*< 


*»»**»*»««» 





FOFETAl 

3 

c* 







• FOFETAl 

4 

c* 

FUNCTION 

FOFETAl COMPUTES THE 

BLADE THICKNESS FOR 

SPECIFIED 

* FOFETAl 

5 

c* 

DISTANCE along THE 

BLADE CHORD 



• FOFETAl 

6 

c» 







• FOFETAl 

7 

c*« 







►** FOFETAl 

8 


DIMENSION 

0ETA(8) , 

COEFFS (51 , 

ETA(2,8) , 


FOFETAl 

9 


1 

X<3», 

xa (3) 

• 



FOFETAl 

10 


Z 

YENO(2> , 

YH(3) 

« 

YSTART(2), 


FOFETAl 

11 


3 

VH (31 





FOFETAl 

12 


CCMMON 

ALFA. 

BETA, 

BLADES, 

8LNTN, 

CBETA, 

FOFETAl 

13 


1 

CH, 

COEFFS, 

DELTAU, 

ETA, 

ETA2, 

FOFETAl 

14 


2 

FTAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

FOFETAl 

15 


3 

IFLAG, 

IREGION, 

J2, 

JFLG, 

KFLG, 

FOFETAl 

16 


4 

NOLAOES, 

OMEGAR, 

OTIME, 

PI, 

R, 

FOFETAl 

17 


5 

rm. 

RMIN, 

S, 

SBETA, 

SNDSPD, 

FOFETAl 

18 


6 

SRCOST, 

SHPDST, 

SHPTAU, 

TAO. 

THKRAT, 

FOFETAl 

19 


7 

TIMTHRU, 

VH, 

X, 

XO, 

YEND, 

FOFETAl 

20 


6 

YH, 

YSTART, 

DELETAl, 

DELF 


FOFETAl 

21 


E = XETA/CH 






FOFETAl 

22 


A = C0F.FFS(1)*SQRT(E) 





FOFETAl 

23 


A1 = C0EFFS(5)*E ♦ COEFFSIAI 




FOFETAl 

24 


Al = Al*E » 

COEFFS(3l 





FOFETAl 

25 


A1 = (A1*E T 

COEFFS (2) ) *E 




FOFETAl 

26 


FOFETAl = (A 

♦ AD’THKRAT 




FOFETAl 

27 


10 return 






FOFETAl 

28 


END 






FOFETAl 

29 


Subroutine OBSPOS 

This subroutine returns X, the three coordinates of the observer relative to. the 
YH-frame fixed to the helicopter. 


Flow Chart of Subroutine OBSPOS 
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Subroutine CALCYST 

Subroutine CALCYST is called by RTN. A second entry point, CALCY, is used by 
subroutine CALCSP. This subprogram calculates the points of intersection of the projec- 
tion of the sphere on the rotor system plane with the circles(s) defined by the rotor sys- 
tem. The calculation is made for either the circle defined by the outer blade radius or 
for, that defined by the hub radius. The two components of the more clockwise point of 
intersection (relative to the observer) are returned in YSTART. The more counterclock- 
wise point is returned in YEND. The angle (at the observer's position) defined by 
YSTART is returned in PHIO. The CALCYST entry contracts the sphere. The CALCY 
entry does not. 
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Listing of Subroutine CALCYST 

SUBROUTINE CALCYST (RSQ, PHIC ) CALCYST 2 

Q, CALCYST 3 

C* • CALCYST 4 

C* SUBROUTINE CALCYST IS CALLED BY PROGRAM RTN. A SECOND ENTRY * CALCYST 5 

C» POINT, CALCY, IS USED BY SUBROUTINE CALCSP. THIS ROUTINE * CALCYST 6 

C* . CALCULATES THE POINTS OF INTERSECTION OF THE PROJECTION * CALCYST 7 

C* OF THE CONTRACTING SPHERE ON THE ROTOR SYSTEM PLANE WITH • CALCYST 0 

C* THE CIRCLE (S) DEFINED BY THE ROTOR SYSTEM. WHEN THE ROUTINE * CALCYST 9 

C* IS ENTERED WITH RSQ EQUAL TO THE OUTER BLADE RADIUS, THE * CALCYST 10 

C* OUTER CIRCLE IS USED. WHEN ENTERED WITH RSQ EQUAL TO THE * CALCYST 11 

C* INNER BLADE RADIUS, THE INNER CIRCLE IS USED. THE X AND Y • CALCYST 12 

C* COMPONENTS OF THE POINT OF INTERSECTION MORE CLOCKWISE * CALCYST 13 

C* (RELATIVE TO THE OBSERVER) ARE RETURNED IN YSTART<2). THE * CALCYST 14 

C* MORE COUNTERCLOCKWISE POINT IS RETURNED IN YENO(2). THE * CALCYST 15 

C* ANGLE (AT THE OBSERVERS POSITION) OEFINEO DY YSTART IS * CALCYST 16 

C» RETURNED IN PHIO. THE CALCYST ENTRY CONTRACTS THE SPHERE. • CALCYST 17 

C* the CALCY ENTRY DOES NOT. * CALCYST 18 

C’ 

C’ 








• CALCYST 

19 







CALCYST 

20 

DIMENSION 

BE TA (8) , 

C0EFFS(5) , 

ETA(2,8), 


CALCYST 

21 

1 

X (3) , 

XO (3) 

t 



CALCYST 

22 

2 

TEND (2) , 

YH(3> 

♦ 

YSTART<2), 


CALCyST 

23 

3 

VH (3) 





CALCYST 

24 

COMMON 

ALFA, 

BETA, 

BLADES, 

BLNTN, 

CBETA, 

CALCYST 

25 

1 

CH, 

COEFFS, 

OELTAU, 

ETA, 

ETA2, 

CALCYST 

26 

2 

ETAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

CALCYST 

27 

3 

IFLAG, 

IREGION, 

J2, 

JFLG, 

KFLG, 

CALCYST 

28 

4 

NBLADES, 

OMEGAR, 

OTIME, 

PI, 

R, 

CALCYST 

29 

5 

RM, 

RMIN, 

S, 

SBETA, 

SNDSPO, 

CALCYST 

30 

6 

SRCDST, 

SWPOST, 

SWPTAU, 

TAU, 

THKRAT, 

CALCYST 

31 

7 

timthru. 

VH, 

X, 

XC, 

YENO, 

CALCYST 

32 

8 

YH, 

YSTAPT, 

OELFTAl, 

OELF 


CALCYST 

33 

IFLAG = 2 






CALCYST 

34 

SHPTAU = SWPTAU ♦ OELTAU 




CALCYST 

35 







CALCYST 

36 

CONTRACT 

THE SPHERE 





CALCYST 

37 







CALCYST 

38 

SWPOST = RM 

- SNDSPD*TIMTHRU»OELTAU 



CALCYST 

39 

IFUFLG.NE. 

0) CALL NEWYH(SWPTAU) 




CALCYST 

40 

CALL OBSPOS 






CALCYST 

41 

OBSANGL = ASIN(X(3)/SQP,T(X(1)**2 

♦ X(2)**2 

♦ X(3)**2)) 


CALCYST 

42 


TEST TERMINATION CRITERIA 

IF (ORSANGL.GT. (PI/4. ) ) GO TO 10 
IF(RM.LT.RMIN) GO TO 10 
ENTRY CALCY 
I FLAG = 1 

SLROSO = X(1)*X(1) + X(2)»X(2) 

OBSOSa = SURDSQ ♦ X(3)*X(3) ' 

DELTA = 0.5*(OBSDSQ *• PSQ - SHPOST*SWPDST ) ' 

DISC = RSQ’SL'ROSO - OELTA*DELTA 
IFtOISC.LE. 3. ) GO TO 10 
IFLAG = 0 

COMPUTE INTERSECTIONS AND INITIAL SHEEP ANGLE 


lu 


FAC = 1.0 
IF(X(1).L 
IF(ABS(X( 
YSTART(2) 
YSTAPT (1) 
YEN0(2) = 
YENO(l) = 
GO TO 7 
YSTART(l) 
YSTART(2) 
YENO(l) = 
YENO(2) = 
SPCOST = 
PHIB = AC 
IF (YSTART 
RETURN 
END 


T.?.) FAC = -l.C 
1)).LT.P) GO TO 5 

= (X(2)»OELTA ♦ X(l)»SQRT(OISC) )/SURDSQ 
= (DELTA - X(2)*YSTART(2) )/X(l) 
(X(2)*OELTA - X (1)*SQRT (DISC) ) /SURDSQ 
(DELTA - X (2) *Y£NO(2) )/X (1) 

= (X(1)*DELTA - FAC»X (2) »DISC) /SURDSQ 
= (OELTA - X(l) »YSTAPT(1) )/X(2) 
(X(l)*OELTA * FAC»X(2)*DISC)/SURDSQ 
(DELTA - X(l) »YENO( 1) ) /X (2) 

SQRT (SWPDST»»2 - X(3)*»2) 

OS(C .9999999999» (YSTART(l) - X ( 1 ) ) /SPCDST ) 
(2).LT.X(2)) PHIO = 2.*PI - PHIO 


CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CAI_CVCT 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 

CALCYST 


43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
6C 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 
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Subroutine CALCSP 

This routine controls the sweeping across the arc defined by a fixed position of the 
contracting sphere. The sweep is in a counterclockwise direction beginning at PHIO and 
terminating when a point is reached outside the circle defined by the rotor system. A 
trapezoidal integration over the intersection of the arc with each blade is performed 
dynamically. The coarseness of the sweep increment is determined as a function of 
whether or not a point is on a blade. Even finer sweep increments are used near the 
leading edge of blades with blunt leading edges. 
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Listing of Subroutine CALCSP 



SUBROUTINE CALCSP (RSQ , POSQ, PHI , Y, 

SCRIPTP) 



CALCSP 

2 

C*' 


«»««««««• 



'• CALCSP 

3 

c* 





• CALCSP 

4 

c» 

SUBROUTINE CALCSP IS CALLED BY 

PROGRAM 

RTN. THIS 

ROUTINE 

* CALCSP 

5 

c* 

CONTROLS the SWEEPING tCOUNTERCLOCKWlSEl ACROSS 

The arc 

* CALCSP 

6 

c* 

OEFINEO BY A FIXED POSITION 

OF THE CONTRACTING 

SPHERE. 

♦ CALCSP 

7 

c* 

TWO sweeping increments are 

USEDI 0ELPHI2 FOR SWEEPING 

* CALCSP 

8 

c* 

across a ROTOR BLADE AND OELPHIl FOR 

THE SPACE 

BETWEEN 

* CALCSP 

B 

c* 

BLADES, AS DETERMTNEO BY SUBROUTINE TEST. THE LINE INTEGRAL 

* CALCSP 

1C 

c» 

REPRESENTING THE CONTRIBUTION OF EACH 

BLADE IS 

A 

• CALCSP 

11 

c* 

TRAPEZOIDAL INTEGRATION. POINTS ON THE EDGE OF 

ANY BLADE 

• CALCSP 

12 

c* 

ARE OBTAINED FROM SUBROUTINES LEFTA AND TEETA. 


• CALCSP 

13 

c* 





* CALCSP 

14 

c*< 





►* CALCSP 

15 


DIMENSION BETA(8>, COEFFS(5), 

ETA(2, 8) , 


CALCSP 

16 


1 X(3>, XQ(3) 


Y<3» , 


CALCSP 

IZ 


2 YEND(21, YH(3l 


YLSTC3), 


CALCSP 

18 


3 YSTARH2), VH(3) 




CALCSP 

19 


COMMON ALFA, BETA, 

BLADES, 

BLNTN, 

CBETA, 

CALCSP 

20 


1 CH, COEFFS, 

OELTAU, 

ETA, 

ETA2, 

CALCSP 

21 


2 ETAMAX, FAC, 

FN, 

FNS, 

FETAl, 

CALCSP 

22 


3 IFLAG, IREGION, 

J2, 

JFLG, 

KFLG, 

CALCSP 

23 


l» NBLAOES, OMEGAR, 

OTIME, 

PI, 

R, 

CALCSP 

24 


5 RM, RMIN, 

S, 

SBETA, 

SNDSPD, 

CALCSP 

25 


6 SRCOST, SWPOST, 

SHPTAU, 

TAU, 

THKRAT, 

CALCSP 

26 


Z TIMTHRU, VH, 


XC, 

YEND, 

CALCSP 

27 


8 YH, YSTART, 

OELETAl, 

OELF 


CALCSP 

28 

c 

■' 




CALCSP 

29 

c 

INITIALIZE VARIABLES FOR THIS 

SWEEP 



CALCSP 

30 

c 





CALCSP 

31 


DELPHII = CH/ (FN*SRCOST> 




CALCSP 

32 


QELPHI2 •= CH/ (FNS»SRCnST) 




CALCSP 

33 


OELETAl = 9. 




CALCSP 

34 


OELF = 0. 




CALCSP 

35 


>LST(3) = C. 




CALCSP 

36 


JREGION = 0 




CALCSP 

37 


SUMSP = 0. 




CALCSP 

38 


OGLST = 0. 



- 

CALCSP 

39 


SCRIOTP = C. 




CALCSP 

40 


J = 1 




CALCSP 

41 


J2 = 1 




CALCSP 

42 

c 





CALCSP 

43 

c 

DETERMINE WHICH REGION HE ARE 

IN 



CALCSP 

44 

c 





CALCSP 

45 


CALL TES T( YSTART, RSQ, ROSQ,SWPTAU» 




CALCSP 

46 


IF (IREGION. EG. 1 > GO TO 20 




CALCSP 

47 


1C IDELREG = IREGION - JREGION 




CALCSP 

48 


IF (lOELREG.EQ.C > GO TO 20 




CALCSP 

49 


IFdDELREG.EO.E) GO TO 3C 




CALCSP 

50 


^ t -r ry r- t r\ t~ r- r> f » ^ f\ -rrs r% n 

J.r»iucur\cA>*c.u«:.z/ w \j lu 





51 


IF (IDELREG. EO. 3) GO TO 130 




CALCSP 

52 


IF(IDELREG.EO.-l)GO TO 80 




CALCSP 

53 

c 





CALCSP 

54 

c 

NOT ON A blade. SELECT LARGE ANGLE AND 

SEARCH The 

ARC 

CALCSP ^ 

55 

c 





CALCSP 

56 


20 DELPHI = DELPHII 




CALCSP 

57 


PHI = PHI ♦ FAC’OELPHI 




CALCSP 

58 


Y(l) = X(l) ♦ SRCPST»COS(PHI( 




CALCSP 

59 


Y(2) = X(2) » SPCDST*SIN(PHI) 




CALCSP 

60 


JREGION = IREGION 




CALCSP 

61 


call TEST(Y,RSQ,R0SQ,SWPTAU) 




CALCSP 

62 


GO TO 19 




CALCSP 

63 

c 





CALCSP 

64 

c 

ENTERING REGION 6 FROM REGION 

G. USE the intersection DATA 

CALCSP 

65 

c 





CALCSP 

66 


3C Y(l) = YSTART(l) 




CALCSP 

67 


Y(2) = YSTART(2( 




CALCSP 

68 

c 





CALCSP 

69 
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ENTERING REGION 6 FROM ANOTHER REGION. SELECT SMALL ANGLE 

CALCSP 

70 


AND PROCESS THE BLADE 

CALCSP 

71 



CALCSP 

72 

4C 

DELPHI = DELPHI2 

CALCSP 

73 


IF(ETA(1,J2).LT.ETAMAX) DELPHI = DELPHI/BLNTN 

CALCSP 

74 


FETAl = FOFETAl (ETA (1, J2) ) 

CALCSP 

75 

50 

JREGION = IREGION 

CALCSP 

76 



CALCSP 

77 


INCREASE THE SWEEP ANGLE AND TEST THE NEXT POINT 

CALCSP 

78 



CALCSP 

79 


PHI = PHI ♦ FAC’DELPHI 

CALCSP 

80 


ETAll = ETA (1,J2) 

CALCSP 

81 


FTA22 = ETA(2,J2) 

CALCSP 

82 


FETAll = EETAl 

CALCSP 

83 


SAVETA = OFLETAl 

CALCSP 

84 


SAVEF = DELE 

CALCSP 

85 


YLST(l) = Yd) 

CALCSP 

86 


YLST(2) = Y(2) 

CALCSP 

87 


Yd) = Xd) ♦ SRCDST*COS(PHI) 

CALCSP 

88 


Y(2) = X(2) SRCDST*SIN(PHI) 

CALCSP 

89 


CALL TEST(Y,RSO,ROSq,SWPTAU) 

CALCSP 

90 


lOELREG = IREGION - JREGION 

CALCSP 

91 


IF(IDELREG.EQ.O) GO TO 6C 

CALCSP 

92 


IF(IOELREG.EO.-5) GO TO 14C 

CALCSP 

93 


IF(IGELREG.EQ.-6) GO TO 170 

CALCSP 

94 


GO TO 180 

CALCSP 

95 



CALCSP 

96 


ON A BLAOE 

CALCSP 

97 



CALCSP 

98 

60 

J = J + 1 

CALCSP 

99 


OELETAl = ETAd,J2) - ETAll 

CALCSP 

100 


0ELETA2 = £TA(2,J2) - ETA22 

CALCSP 

ICl 


FETAl = EOFETAl (ETA (1, J2) ) 

CALCSP 

102 


DELE = FETAl - EETAH 

CALCSP 

103 


OELGAM = SQRT(0ELF**2 DELETA1»»2 ♦ DELETA2**2) 

CALCSP 

1C4 


IF(0ELETA1.LT.C.) OELGAM = - OELGAM 

CALCSP 

1C5 


SP = PFUNCN(£TA22,OELF, OELETAl, YLST) 

CALCSP 

1C6 


SUMSP = SUMSP ♦ SP*0ELGAM 

CALCSP 

107 


SP = PFUNCN (ETA22, SAVEF, SAVETA, YLST) 

CALCSP 

ICB 


SUMSP = SUMSP ♦ SP’OGLST 

CALCSP 

109 


OGLST = OELGAM 

CALCSP 

110 


1F(IREGICN.EQ,6) GO TO 50 

CALCSP 

111 



CALCSP 

112 


JUST LEFT REGION 6. PERFORM THE TFAPEZOIDAL INTEGRATION 

CALCSP 

113 



CALCSP 

114 


SP = PFUNCN (ETA (2, J2) ,OELF, OELETAl, Y) 

CALCSP 

115 


SUMSP = SUMSP * SP»0ELGAM 

CALCSP 

116 


IFdPEGION.EQ.C) GO TO 80 

CALCSP 

117 


OGLST = 0 . 

CALCSP 

. 118 


DELETAl = C. 

CALCSP 

119 


DELE = C. 

CALCSP 

120 


J = 1 

CALCSP 

121 


J2 = 1 

CALCSP 

122 


GO TO 20 

CALCSP 

123 



CALCSP 

124 


PASSED OUT OF THE ROTOR SYSTEM FOP THIS ARC 

CALCSP 

125 



CALCSP 

126 

80 

SCRIPTP = 0.5*nELTAU*SUMSP/SHPDST 

CALCSP 

127 


RE TURN 

CALCSP 

128 



CALCSP 

129 


ENTERING REGION 6 FROM REGION 1 

CALCSP 

130 



CALCSP 

131 

90 

IF(ETA (1,J2) .GT.CH/2.) GO TO 100 

CALCSP 

132 



CALCSP 

133 


THROUGH THE LEADING EDGE 

CALCSP 

134 



CALCSP 

135 


CALL LEETA(ETA(2,J2),YI 

CALCSP 

136 


FTA(1,J2) = 0. 

CALCSP 

137 


GO TO 110 

CALCSP 

138 



CALCSP 

139 


THROUGH THE TRAILING EDGE 

CALCSP 

14C 


40 
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100 CALL TEETA(ETA(Z,J2»,Y) 

ETA« 1, J2) = CH . 

110 ETA(2,J2) = ETA2 

SET ANCLE AND RETURN TO BLADE LOGIC 

120 PHI = AC0S(J.B99gg99990»(Y(l) - X ( 1 » » /SRCDST) 

IF(Y(2).LT.X(2) ) PHI = 2.*PI - PHI 
GO TO 40 

ENTERING REGION 6 FROM REGION 3. START FROM INTERSECTION 
WITH INNER RADIUS 

130 CALL CALCY (ROSQ.PHX) 

Y(l) = YENDCl) 

Y(2I = YEN0C2) 

ETA(1,J2) = CH/2. - Y{1I»SBETA ♦ Y(2)»CBETA 

IF (ETA (1, J2) .LT. C. ) GO TO 90 

ETA(2,J2) = -Y(1)*CBETA - Y(2)*SBETA 

GO TO 120 - , 

ENTERING REGION 1 FROM REGION 6 

140 IF (ETA (1, J2) .GT.CH/2. ) GO TO 150 

THROUGH LEADING EDGE 

call LFETA(ETA(2,J2I.Y) '' ' ' 

FTA(1,J2) = C. 

GO TO 16C ' ■ ' 

THROUGH TRAILING EDGE 

150 call TEETA(ETA(2,J2»,Y) 

ETA(1,J2) = CH 
16C ETA(2,J2I = ETA2 
G 0 T 0 6 0 

LEAVING the ROTOR SYSTEM FOR THIS ARC. GET INTERSECTION WITH 
OUTER RADIUS 

170 CALL CALCY (RSQ.PHXI 
Y (1) = YENDdl 
Y(2) = YENO(2) 

GO TO 190 

ENTERING REGION 3 FROM REGION 6. GET ' INTERSECT ION WITH 
INNER RADIUS 

180 CALL CALCY(ROSQ.PHX) 

Y(l) = YSTART(I) 

Y(2) = YSTARKF) 

190 ETA(1,J2) = CH/2. - Y(1)*SBETA ♦ Y(2I»CBETA 
IF (ETA (1, J2I .LTi 0. I GO TO 140 
ETA(2,J2) =■ -Y(1(»CBETA - Y(2I*SBETA 
GO TO 6C ■ 

END • ' ■ 


CALCSP 

141 

CALCSP 

142 

CALCSP 

143 

CALCSP 

144 

CALCSP 

145 

CALCSP 

146 

CALCSP 

147 

CALCSP 

148 

CALCSP 

149 

CALCSP 

15C 

CALCSP 

151 

CALCSP 

152 

CALCSP 

153 

CALCSP 

154 

CALCSP 

155 

CALCSP 

156 

CALCSP 

157 

CALCSP 

158 

CALCSP 

159 

CALCSP 

160 

CALCSP 

161 

CALCSP 

16 2 

CALCSP 

163 

CALCSP 

164 

CALCSP 

165 

CALCSP 

166 

CALCSP 

167 

CALCSP 

16 8 

CALCSP 

169 

CALCSP 

17C 

CALCSP 

171 

CALCSP 

172 

CALCSP 

173 

CALCSP 

174 

CALCSP 

175 

CALCSP 

176 

CALCSP 

177 

CALCSP 

178 

CALCSP 

179 

CALCSP 

180 

CALCSP 

181 

CALCSP 

182 

CALCSP 

183 

CALCSP 

184 

CALCSP 

185 

CALCSP 

186 

CALCSP 

187 

CALCSP 

188 

CALCSP 

189 

CALCSP 

190 

CALCSP 

191 

CALCSP 

192 

CALCSP 

193 

CALCSP 

194 

CALCSP 

195 

CALCSP 

196 

CALCSP 

197 

CALCSP 

198 
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Function PFUNCN 


This function subprogram computes the integrand of the line integral at each sweep 
position. 


Flow Chart of Function PFUNCN 



Listing of Function PFUNCN 



FUNCTION PFUNCN(XETA,0 

elf.deletai. 

Y) 



PFUNCN 

2 



PFUNCN 

3 

c» 





* 

PFUNCN 

4 

c* 

FUNCTION PFUNCN IS 

called by subroutine 

CALCSP. THl 

:$ FUNCTION • 

PFUNCN 

5 

c* 

COMPUTES THE INTEGRAND OF THE 

LINE integral at EACH * 

PFUNCN 

6 

c* 

SWEEP POSITION. 




* 

PFUNCN 

7 - 

c* 





* 

PFUNCN 

8 


PFUNCN 

g 


DIMENSION BETA (S» , 

C0EFFS(5» , 

ETA(2,8> , 


PFUNCN 

10 


1 X (3> , 

X0(3», 




PFUNCN 

11 


2 Y£N0(2), 

YH(3» , 


YSTART(2», 


PFUNCN 

12 


-3 VH<3) 





PFUNCN 

13 


DIMENSION Y(3) 





PFUNCN 

14 


COMMON ALFA, 

BETA, 

BLADES, 

BLNTN, 

CBETA, 

PFUNCN 

15 


1 CH, 

COEFFS, 

OELTAU, 

ETA, 

ETA2, 

PFUNCN 

16 


2 ETAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

“FUNCN 

17 


3 IFLAG, 

IREGION, 

J2, 

JFLG, 

KFLG, 

PFUNCN 

18 


4 NBLADES, 

OMEGAR, 

OTIME, 

PI, 

R, 

PFUNCN 

IB 


5 PM, 

RMIN, 

S, 

SBETA, 

SNDSPD, 

PFUNCN 

20 


6 SPCOST, 

SHPOST, 

SHPTAU, 

TAU, 

THKPAT, 

PFUNCN 

21 


7 TIMTHRU, 

VH, 

X, 

XC, 

YENO, 

PFUNCN 

22 


8 yh,ystart,oeletax,oelfx 




PFUNCN 

23 


PFUNCN = 0. 





PFUNCN 

24 


S = (((X(l) - Y(U)*SBETA ♦ (Y(2) 

- X(2K» 

CBETA)*DELF 


PFUNCN 

25 


1 ♦ (X(3> - Y(3n*OELETAl)/SHPOSr 




PFUNCN 

26 


Q = OELF**2 t DELETA1**2 - S»*2 




PFUNCN 

27 


IFCQ.LT. I.E-IOI RETURN 





PFUNCN 

28 


PFUNCN = OELF* (OMEGAR»XETA - VH(2» 

*CBETA ♦ 

VH(1»*SBETA>/SQRT(Q) 

PFUNCN 

2B 


RETURN 





PFUNCN 

30 


END 





PFUNCN 

31 


Subroutine LEETA 

Subroutine LEETA returns the coordinates (in the YH frame) of the point of inter- 
section of a fixed arc and the leading edge of a blade. The second ETA component (see 
fig. 2) is also returned. A second entry point, TEETA, does the same thing for a trailing 
edge. 
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Flow Chart of Subroutine LEETA 
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Listing of Subroutine LEETA 



SUBROUTINE LEETA(XETA, 

Y( 




LEETA 

2 


3 

c* 







* LEETA 

4 

c* 

SUBROUTINE LEETA IS 

called by subroutine 

CALCSP. 

THIS 

* LEETA 

5 

r* 

ROUTINE 

RETURNS THE COORDINATES (IN THE Y REFERENCE 

* LEETA 

6 

c* 

FRAME) 

OF THE INTERSECTION OF 

A FIXED 

SWEEPING 

ARC AND 

* LEETA 

7 

c* 

The leading edge 

OF A blade. 

A SECOND 

ENTRY POINT, TEETA, 

* LEETA 

8 

c» 

DOES THE SAME FOR 

THE TRAILING EDGE. THE SECDND 

ETA 

* LEETA 

9 

c* 

COMPONENT IS ALSO 

RETURNED. 




• LEETA 

10 

c» 







* LEETA 

11 

c»»* 

****»♦»*»»»»»»*♦»*»****»♦•**»»*•♦»•*•*****»*********•*♦*•*********** LEETA 

12 


DIMENSION 

BETA (8) , 

COEFFS (5) . 

ETA (2, 8) , 


LEETA 

13 


1 

X(3) , 

XO (3) , 


Y ( 3 ) , 


LEETA 

14 


Z 

YEND(2), 

YH(3( , 


YSTART (2) 


LEETA 

15 


3 

VH (3) 



■ ..... 


LEETA 

16 


COMMON 

ALFA, 

BETA, 

BLADES, 

BLNTN, 

CBETA, 

LEETA 

17 


1 

CH, 

COEFFS, 

OELTAU, 

ETA, 

ETA2, 

LEETA 

18 


2 

ETAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

LEETA 

19 


3 

IFLAG, 

IREGION, 

J2, 

JFLG, 

KFLG, 

LEETA 

2C 


4 

NBLAOES, 

OMEGAR, 

OTIME, 

PI, 

P, 

LEETA 

21 


5 

RM, 

RMIN, 

S, 

SBETA, 

SNDSPD, 

LEETA 

22 


6 

SRCOST, 

SWPOST, 

SHPTAU, 

TAU, 

THKRAT, 

LEETA 

23 


7 

TIMTHRU, 

VH, 

X, 

xc,.. 

YENO, 

LEETA 

24 


fl 

YH, 

YSTA.RT, 

OELETAl, 

DELF 


LEETA 

25 


A = -0.5*CM 






LEETA 

26 


GO TO 1C 






LEETA 

27 


ENTRY TEETA 






LEETA 

28 


A = 0.5*CH 






LEETA 

29 

1 

0 R = -A 






LEETA 

30 


XBARl = X(l) 

♦ A'SBETA 





LEETA 

31 


XBAR2 = X(2) 

♦ B»C0ETA 





LEETA 

32 


XBAR = X(l)» 

CBETA ♦ X( 

2)*SBETA 




LEETA 

33 


XTILOE = XBAR1*XRAP1 + 

XBAR2*XOAR2 

- SRCDST 

•SRCDST 


LEETA 

34 


DISCI = SQRT(X0AR*XBAR 

- XTILOE) 




LEETA 

35 


ETA2 = -XBAR 

♦ DISCI 





LEETA 

36 


IF (ABS (XETA 

- ETA2) ,GE 

.CH/FN) ETA2 

= ETA2 - 

2,*DISC1 


LEETA 

37 


Yfl) = a»SBETA -ETA2»CBETA 




LEETA 

38 


Y(2) = A»CBETA -FTA2*S8ETA 




LEETA 

39 


RETURN 






LEETA 

40 


END- 






LEETA 

41 


Subroutine TEST 

TEST returns the region index (IREGION) for a given pair of coordinates in the YH 
reference frame at a given time. The region indices correspond to points outside the 
outer radius, inside the hub radius, within the rotor system but not a blade, and on a blade. 
When a point is found to be on a blade, the routine returns the blade number, the angle 
made by that blade, the sine and cosine of this angle, and the ETA coordinates of the point. 
The subroutine included with the current version of RTN is limited to rotor systems with 
rectangular blades. 
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Listing of Subroutine TEST 


c*** 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c» 

c» 

c* 

c* 

c» 

c* 

c* 

c» 

c* 

c* 


SUBROUTINE TEST ( Y , RSQ , ROSQ, TI HE ) 

SUBROUTINE TEST IS CALLEO BY SUBROUTINE CALCSP. THIS 
ROUTINE RETURNS THE REGION INDEX EOR A GIVEN TIME AND 
A GIVEN PAIR OF COORDINATES IN THE Y REFERENCE FRAME. 
THE REGION INDEX HAS FOUR POSSIBLE VALUES. 


IREGION = Cl 


= It 


THE POINT lies OUTSIDE THE CIRCLE 
DEFINED BY THE BLADES OUTER RADIUS 

THE POINT LIES BETWEEN THE CIRCLES 
DEFINED BY THE BLADES OUTER AND 
INNER RADII, BUT IS NOT ON A BLADE 

THE POINT LIES INSIDE THE CIRCLE 
DEFINED BY THE BLADES INNER RADIUS 


= 61 THE POINT LIES ON A BLADE., 

WHEN THE REGION INDEX IS 6, TEST RETURNS THE BLADE NUMBER, 
J2, THE ANGLE MADE BY THAT BLADE, THE SIN AND COSINE OF 
THIS ANGLE, SBETA AND CBETA, AND THE ETA COORDINATES OF 
THE POINT. 


DIMENSION 


COMMON 




TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST • 

TEST 

TEST 

TEST 


RYSa = Y (II **2 ♦ Y(2) *»2 
IREGION = C 

lE(RYSQ.GT.RSQ) GO TO 3 
IREGION = 3 

IF(RYSQ.LT.ROSQI GO TO 3 
IREGION = 1 

LIMTOP = J2 ♦ NBLAOES - 1 
DO 1 J=J2, LIMTOP 
I = J 

IF (J.GT.NBLADES) I = J 
BETA(I) = OMEGAR*TIME 
CPETA = COS(BETA(I) I 
SBETA = SIN(BETA(I) ) 

IF((Y(2)*CBETA - Y ( 1 1 ♦SBE T A ) •• 2 . GT . CH’CH/' 
ETA(2,I) = -Y (1 >»CBETA' - Y(21*S8ETA 
IF (ETA (2, 1 1 ,GT. C . 1 GO TO 2 
CONTINUE 

COS(BETA(J211 
SIN(BETA(J2)I 


- NBLAOES 

+ FLOAT (I-l) *6. 2fl31 


8E 30 8/FLOAT (NBLAOES) 
L. ) GO TO 1 


CBETA = 
SBETA = 
GO TO 3 
IREGION 
J2 = I 
ETAd, J2) 
RETURN 
END 


= 6 


= CH/2. - Y(1)*SBETA ♦ Y(2I*C8ETA 


TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST . 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 

TEST 


2 

3 

<» 

5 

6 

7 

8 
9 

1C 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


BETA(8), 

COEFFS (5) , 

ETA(2,8) , 


TEST 

28 

X(3), 

XO (3) 

t 

Y(3) , 


TEST 

29 

YEND(2) , 

YH(3) 

♦ 

YSTART(2) , 


TEST 

30 

VH(3) 





TEST 

31 

ALFA, 

BETA, 

BLADES, 

BLNTN, 

CBETA, 

TEST 

32 

CH, 

COFFFS, 

OELTAU, 

ETA, 

ETA2, 

TEST 

33 

ETAMAX, 

FAC, 

FN, 

FNS, 

FETAl, 

■ TEST 

34 

IFLAG,- 

IREGION, 

J2. 

JFLG, 

KFLG, 

TEST 

3.5 

NBLAOES, 

OMEGAR, 

OTIHE. 

PI, 

P, 

TEST 

36 

RM, 

RMIN, 

s, ■ 

SBETA, 

SNOSPD, 

TEST 

37 

SRCDST, 

SHPOST, 

SHPTAU, 

TAU, 

THKRAT, 

• TEST 

38 

TIMTHRU, 

VH, 

X, , 

XD, 

YEND, • 

■ TEST. 

39 

YH, 

YSTART, 

OELETAl, 

OELF 


. TEST 

40 


41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 
'5 3 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 
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APPENDIX B 


ROTOR THICKNESS NOISE POSTPROCESSOR PROGRAM 

This appendix includes brief descriptions of each of the RTNPP subprograms. 
Where possible, each description is followed by a flow chart and FORTRAN listing of the 
subprogram. The CSDS subroutine is the property of International Mathematical and 
Statistical Library (IMSL). Four other subroutines, SPLDER, PSEUDO, INFOPLT, and 
CALPLT, are the property of Langley Research Center. For these five subroutines, 
FORTRAN listings are not included, but the more detailed descriptions (ref. 7) provided 
should permit their replacement by similar routines which might be available at other 
facilities. 


Program RTNPP 

Program RTNPP reads the output from the RTN program. If the contents of 
TAPE 12 were disposed to the card punch following execution of the RTN program, this 
deck should be placed in the RTNPP input stream. If TAPE12 is a disk file, it should be 
rewound prior to executing RTNPP. No other input is necessary for the postprocessor 
program. Once the data has been read and stored into arrays, trigonometric smoothing 
is used to obtain an estimate of the validity of each data point in the T>(X,t) plotted against 
t spectrum. These estimates are then passed to a cubic spline smoothing routine. Inter- 
polation is performed on the smoothed data to provide values at equally spaced values of 
the observer time. The smoothed data are numerically differentiated with respect to the 
observer time and the results are plotted. 
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Listing of Program RTNPP 



PROGRAM RTNPP«INPUT=101B,OUTPUT=1C1B, TAPE5= INPUT, TAPE 6=0 UTPUT ) 

RTNPP 

2 


DIMENSION SCRIPTP (35C ) , OT IME ( 35 0) , F (500 » , T < SCO ) 

,OER1(500) , 

RTNPP 

3 


1 CER2 ( 500 1 ,COEF (350,4) ,HK (38501 , OF (350 


RTNPP 

4 


OIMENSION XO (3) ,COEFFS(5l 


RTNPP 

5 


LOGICAL MOROATA 


RTNPP 

6 


NAMELIST/ROTOR/XC ,NBLADES,R, CH, OMEGA, RO.THKR AT, 

OELTET,N,NS,KFLG, 

RTNPP 

7 


1 ETAMAX,BLNTN,JFLG,COEFFS,TAUINT,PEROM,OTAUM,SNDSPD, MOROATA 

RTNPP 

8 

1 

FORMAT (1hC,*IERR = *,I5) 


RTNOP 

9 

2 

FORMAT (15, 3E25. 15 ) 


RTNPP 

10 


M i 0 


RTNPP 

11 


INIT = 0 


RTNPP 

12 

1C 

PEAD(5, ROTOR) 


RTNPP 

13 


HRITE(6, ROTOR) 


RTNPP 

14 

20 

CONTINUE 


RTNPP 

15 


M = M ♦ 1 


RTNPP 

16 

30 

CONTINUE 


RTNPP 

17 


READ (5,2) I,7,0TIME(M),SCR1PTP(M) 


RTNPP 

18 


OTIME(M) = ICCG .*0TIME(M) 


RTNPP 

19 


IF(I.NE.O) GO TO 2C 


RTNPP 

20 


M = M - 1 


RTNPP 

21 


MNPTS = 35C 


RTNPP 

22 


NCVS = 1 


RTNPP 

23 


M = M/3 


RTNPP 

24 


M = 3»M 


RTNPP 

25 


OT = (OTIME(M( - OTIME(I) )/499. 


RTNPP 

26 


00 40 IJ=1,5C0 


RTNPP 

27 


T(IJ) = OTIME(l) ♦ FLOAT(IJ-l) *DT 


RTNPP 

28 

<»0 

CONTINUE 


RTNPP 

29 


MM = 500 


RTNPP 

30 


S = FLOAT(M) ♦ SORT (FLOAT (2*M ) ) 


RTNPP 

31 


PI = 3.141592654 


RTNPP 

32 


NMl = M - 1 


RTNPP 

33 


00 50 1=1, M 


RTNPP 

34 


OERl(I) = 0. 


RTNPP 

35 


HK(I) = 0. 


RTNPP 

36 


DFd) = SCRIPTPd) - SCRIPTP(l) - FLO AT ( I- 1 ) • { SCR IPTP (M) - SCRIPTP 

RTNPP 

37 


1 (1) ) /FLOAT (NMl ) 


RTNPP 

38 

50 

CONTINUE 


RTNPP 

39 


NNMM= MINC(85,M) 


RTNPP 

40 


no 70 J=1,NNMM 


RTNPP 

41 


OC 60 1=2, NMl 


RTNPP 

42 


WK(J) = HK(J) ♦ DF( I) *SIN (FLOAT (J* (I-l) ) *PI/FLOAT (NMl ) ) 

PTNPP 

43 

6C 

CONTINUE 


RTNPP 

44 


WKU) = 2.*HK(J)/FLCAT(NM1) 


RTNPP 

45 

7 0 

CONTINUE 


RTNPP 

46 


DO 9C 1=2, NMl 


RTNPP 

47 


PC so J=1,NNMM 


RTNPP 

48 


OERKI) = OEPKI) ♦ HK(J) *SIN (FLOAT (J*(I-1) )*PI/FLOAT (NMl) ) 

RTNPP 

49 

SC 

CONTINUE 


RTNPP 

50 


COEF(I,l) = OEPKI) ♦ SCRIPTPd) ♦ FLOAT(I-l)*( 

SCRIPTP (M) 

RTNPP 

51 


1 - SCRIPTPd) )/FLOAT(NMl) 


RTNPP 

52 

90 

CONTINUE 


RTNPP 

53 


CCEF(1,1) = SCRIPTPd) 


RTNPP 

54 


COEF(M,l) = SCRIPTP(M) 


RTNPP 

55 


00 100 1=1, M 


RTNPP 

56 


DFd) =1.5*ABS(COEF(I,l) - SCRIPTP(I)) 


RTNPP 

57 

100 

CONTINUE 


RTNPP 

58 


ML = 0 


RTNPP 

59 


IFdNIT.GT.D) ML = 1 


RTNPP 

60 


INIT = INIT ♦ 1 


RTNPP 

61 


CALL PLOTPT (OTIME,SCRIPTP,M,ML, 1) 


RTNPP 

62 


ML = 1 


RTNPP 

63 


IW = -1 


RTNPP 

64 
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CALL CSOS (MNPTS,M,OTIME,SCRIPTP, OF,S, IH, COEF.WK, IERR> 

RTNPP 

65 

WRrTE(6.1» lERR 

RTNPP 

66 

CALL PLOTPT tOTIME.COEF.M.ML,?) 

RTNPP 

67 

XMN = OTIME (U 

RTNPP 

68 

XMX = OTIME (M) 

RTNPP 

69 

YMN = I.EIO 

RTNPP 

70 

YMX = -l.ElO 

RTNPP 

71 

00 110 1=1. M 

RTNPP 

72 

IF(SCRIPTP(I) .LT.YMNI YMN = SCRIPTPdt 

RTNPP 

73 

IFfSCRIPTPai .GT. YMXI YMX = SCRIPTPfl) 

RTNPP 

74 

IF (COEF ( I, 1) .LT.YMNI YMN = C0EF<I,1I 

RTNPP 

75 

IF(C0EF(I,1).GT.YMXI YMX = C0EF<I,1) 

RTNPP 

76 

110 continue 

RTNPP 

77 

0TIME(M+1I = XMN 

RTNPP 

78 

0TIME(M*2I = XMX 

RTNPP 

79 

SCRIPTP<M+1) = YMN 

RTNPP 

80 

SCRIPTP(M*2I = YMX 

RTNPP 

81 

COEF(M+l,ll = YMN 

RTNPP 

82 

COEF(M+2,l) = YMX 

RTNPP 

83 

CALL PLOTPT (OTIME , SCP IPTP ,M .ML , 41 

RTNPP 

84 

CALL PLOTPT (OTIME. COEF, M. ML. 51 

RTNPP 

85 

IN = -1 

RTNPP 

86 

CALL SPL0ER(MNPTS,M,NCVS,MM,MM,0TIME,C0EF,T,F,0ER1.DER2,IH,HK, 

RTNPP 

87 

1 lERR) 

RTNPP 

88 

HRITE(6,1) lERR 

RTNPP 

89 

00 120 1=1, MM 

RTNPP 

90 

OERKII = lOCO.’OERKII 

RTNPP 

91 

120 CONTINUE 

RTNPP 

92 

ML = 2 

RTNPP 

93 

IF(MOROATA) ML = 1 

RTNPP 

94 

CALL PLOTPT(T,OERl,MM,ML,3) 

RTNPP 

95 

M = 0 

RTNPP 

96 

IF(MOROATA) GO TO 10 

RTNPP 

97 

STOP 

RTNPP 

98 

END 

RTNPP 

99 


Subroutine PLOTPT 

This subroutine serves as an interface between RTNPP and the graphics output 
library. Plot size and captions are provided for each plot. 
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Listing of Subroutine PLOTPT 



SUBROUTINE PLOTPT (X,Y,N, ML. J) 


PLOTPT 

2 

DIMENSION XClI.Yd) 


PLOTPT 

3 

IF(ML.EQ.0» CALL PSEUDO 


PLOTPT 

4 

IF(J.GT.3I GO TO 15 


PLOTPT 

5 

XMIN = I.EIQ 


PLOTPT 

6 

YHIN = I.EIO 


PLOTPT 

7 

XMAX =-l.E10 


PLOTPT 

8 

YMAX =-l.E10 


PLOTPT 

9 

00 10 I=1,N 


PLOTPT 

10 

IF (X (I) .LT. XMIN) XMIN = X(I) 


PLOTPT 

11 

IF(Y(I).LT.YMIN) YMIN = Y(I) 


PLOTPT 

12 

IF (X «I) . GT. XMAX) XMAX = X(I) 


PLOTPT 

13 

IF(Y(I).GT.YMAX) YMAX = Y<I) 


PLOTPT 

14 

10 CONTINUE 


PLOTPT 

15 

GO TO 16 


PLOTPT 

16 

15 CONTINUE 


PLOTPT 

17 

XMIN = X(N+1) 


PLOTPT 

18 

XMAX = XlN+2) 


PLOTPT 

19 

YMIN = Y(N+1) 


PLOTPT 

20 

YMAX = YIN»2) 


PLOTPT 

21 

16 CONTINUE 


PLOTPT 

22 

SH = 8. 


PLOTPT 

23 

SW = 5. 


PLOTPT 

24 

ISH = 8 


PLOTPT 

25 

ISV = 5 


PLOTPT 

26 

CALL HELP(XMIN, XM AX , I SH , XMI N . XMAX ) 


PLOTPT 

27 

CALL HELP (YMIN, YMAX, ISV, YMIN, YMAX) 


PLOTPT 

28 

GO TO (20, 30, 40, 60, 70), J 


PLOTPT 

29 

20 CALL INF0PLT(1,N,X,1,Y,1, XMIN, XMAX, YMIN, YMAX, 1., 


PLOTPT 

30 

1 13, 13HOBSERVEP TIME, 18, 18HSCRIPTP (RAH OATA),C,SH, 

SV, l.,l.) 

PLOTPT 

31 

GO TO 50 


PLOTPT 

32 

30 CALL INFOPLTd, N,X,1,Y,1, XMIN, XMAX, YMIN, YMAX, 1., 


PLOTPT 

33 

1 13, 13H0BSERVER T IME , 16 , 1 6HSCRIPTP (SPL INE) , 0 , SH , SV 

, 1 • , 1 , ) 

PLOTPT 

34 

GO TO 50 


PLOTPT 

35 

40 CALL INFOPLTd, N,X,l,Y,t, XMIN, XMAX, YMIN, YMAX, 1., 


PLOTPT 

36 

1 28.28H0BSERVER TIME ( M ULI SECONDS) ,1 7, 1 7HPRESSURE 

(N/M**2) ,0 , 

PLOTPT 

37 

2 SH,SV,1,,1.) 


PLOTPT 

38 

GO TO 50 


PLOTPT 

39 

60 CALL INFOPLT(0,N,X,l,Y,l, XMIN, XMAX, YMIN, YMAX, 


PLOTPT 

4C 

1 1., 13, 13HOBSERVER T I ME , 7 , 7HSCR IPTP , 0 , SH, SV , 1 . , 1 . ) 


PLOTPT 

41 

GO TO 50 


PLOTPT 

42 

70 CALL INFOPLTd, N,X,1,Y,1, XMIN, XMAX, YMIN, YMAX, 


PLOTPT 

43 

1 1., 13, 13H0BSERVER TI ME ,7 , 7HSCR IPTP, 0 , SH , SV , 1 . , 1 . ) 


PLOTPT 

44 

50 IF(ML.EQ.2) CALL CALPLT (0 . , 0 . , 999) 


PLOTPT 

45 

RETURN 


PLOTPT 

46 

END 


PLOTPT 

47 


Subroutine HELP 

This routine serves to eliminate the problem of rounding associated with the anno- 
tation of the values of the variable at the tic marks. If the specified values for the mini- 
mum and maximum result in rounding, a smaller minimum and larger maximum are sub- 
stituted so that rounding does not occur. 
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Flow Chart of Subroutine HELP 



Listing of Subroutine HELP 


SUBROUTINE HELP (XI, X 2.0, XNIN, XM4XI 
INTEGER 0 

PI = ALOGlO <ARS (XI) ) 

R2 • AL OG 1 G ( H S3 ? X t ) J 
P = AMAXl (P1.P2) 

IF(P.LT.O.) P = P - 1. 

IP = P - 2. 

XIP = X1/1C.**IP 
X2P = X2/10.**IP 


IX IP 

= 

XIP - 1 

• 



IX2P 

= 

X2P ♦ 1 

• 



IDIF 

= 

IX2P - 

IXIP 



N = 

lOI 

F/ (0*1) 




JDIF 

= 

N* (O+l) 




IF(I 

OIF 

•EQ.JOIF) GO 

TO 

10 

N = 

N + 

1 




KDIF 

= 

N*(0*1I 




LOIF 


(KOIF - 

lOIF 

)/2 


IXIP 

= 

IXIP - . 

LOIF 



1X2P 

= 

IXIP ♦ 

KOIF 



ID CONT 

INUE 




XMIM 

= 

FLOATCI 

XIP)* 

10 . 

**IP 

XMAX 

= 

FLOATd 

X2P)* 

10. 

**IP 


RETURN 

END 


HELP 

2 

HELP 

3 

HELP 

4 

HELP 

5 

HELP 

6 

HEL° 

7 

HELP 

8 

HELP 

9 

HELP 

10 

HELP 

11 

HELP 

12 

HELP 

13 

HELP 

14 

HELP 

15 

HELP 

16 

help 

17 

HELP 

18 

HELP 

19 

HELP 

2C 

HELP 

21 

HELP 

22 

HELP 

23 

HELP 

24 

HELP 

25 

HELP 

26 
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Subroutine CSDS 

Language : FORTRAN 

Purpose : To fit a smooth cubic spline curve to a set of data points representing a 

univariate function. The functional values are allowed to vary by specified amounts in 
smoothing the curve. Data points may be unequally spaced. 

Use : CALL CSDS (MAX,IX,X,F,DF,S,IPT,COEF,WK,IERR) 


MAX 

IX 

X 

F 

DF 

S 

IPT 

COEF 


An input integer specifying the maximum number of data points for the 
independent variable, as given in the dimension statement of the calling 
program. 

An input integer specifying the actual number of data points for the inde- 
pendent variable. IX = MAX. 

A one -dimensional input array dimensioned at least IX in the calling program. 
Upon entry to CSDS, X(I) must contain the I-th value of the independent 
variable. 

A one-dimensional input array dimensioned at least IX in the calling program. 
Upon entry to CSDS, F(I) must contain the I-th value of the function. 

A one -dimensional input array dimensioned at least IX in the calling program. 
Upon entry to CSDS, DF(I) must contain an estimate of the standard deviation 
of F(I). The value DF(I) determines the amount of variation permitted in the 
I-th functional value in smoothing the data. If the standard deviations of the 
functional values cannot be estimated, specifying DF = 1 is suggested. 

A nonnegative input parameter which controls the extent of smoothing. A 
value in the rar^e (IX - (2*IX)**.5) = S = (IX + (2*IX)**.5) is suggested. 
Increasing values of S permit increasing amounts of smoothing of the func- 
tional values. If S = 0 is specified, an unsmoothed cubic spline curve pass- 
ing through the data points is computed. 

An input initialization parameter. The user must specify IPT = -1 when- 
ever a new X-array is input. The routine will then check to insure that the 
X-array is in strictly increasing order. 

A two-dimensional output array which must be dimensioned (MAX,4) in the 
calling program. Upon return, COEF(I,J) contains the J-th coefficient of the 
spline for the interval beginning at point X(I). The functional value of the 
spline at abscissa XI, where X(I) = XI = X(I + 1) is given by: 
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F(X1) = ((C0EF(I,4)*H + C0EF(I,3))*H + COEF(I,2))*H + C0EF(I,1) 
where H = XI - X(I) 

WK A one-dimensional work area array dimensioned at least (7*IX + 9) in the 
calling program. 

lERR An output integer error parameter. 

= 0 Normal return. 

= J The J-th element of the X-array is not in strictly increasing order. No 
calculation performed. 

= -1 There are less than four values in the X-array. No calculation 
performed. 

Upon return, this parameter should be tested in the calling program. 


Method : A set of IX data points (X^.F^), i = 1,2,. . .IX are given, where F^ is the 

functional value at point Xj. A set of weights DFi,DF2. . -DFix and a nonnegative 
smoothing parameter S are also given. 

The routine computes a cubic spline G(X) having the following properties: 




dx = Minimum of all splines satisfying property 


(i) 


Subroutine SPLDER 

Language : FORTRAN 

Purpose : To perform a cubic spline approximation ; interpolation, and differentiation. 

SPLDER computes F = fi(X), DERI = f[(X), and DER2 = f^’(X) for any number of 
different dependent variable arrays associated with the independent variable array. 

Use : CALL SPLDER (MNPTS,N,NCVS,MMAX,M,X,Y,T, F, DER1,DER2,IW,WK, 

lERR) 

MNPTS An input integer specifying the maximum number of values in the independent 
variable array as stated in the dimension statement of the calling program. 
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N An input integer specifying the number of values in the independent variable 

array. N = MNPTS. 

NCVS An input integer specifying the number of dependent variable tables associated 
with the independent variable. 

NMAX An input integer specifying the maximum number of values at which interpola- 
tion is desired as stated in the dimension statement of the calling program. 

M An input integer specifying the number of values at which interpolation is 

desired on this entry into SPLDER. M = MMAX. 

X A one-dimensional input array containing the independent variables. The 

array X should be dimensioned by at least N in the calling program, and the 
value must be monotonic. 

Y A two-dimensional input array containing the dependent variables. The array 

Y is dimensioned with variable dimension in the subroutine; therefore, Y must 
be dimensioned in the calling program with first dimension MNPTS and sec- 
ond dimension at least NCVS. 

T A one-dimensional input array containing the values of the independent vari- 

able for which values of the dependent variable and the first and second 
derivatives are desired. The array T must be dimensioned by at least M in 
the calling program. 

F A two-dimensional output array in which SPLDER stores the values of the 

function at the M values of the independent variable. The array F is dimen- 
sioned with variable dimension in the subroutine; therefore F must be dimen- 
sioned in the calling program with first dimension MMAX and second dimen- 
sion at least NCVS. 

DERI A two-dimensional output array in which SPLDER stores the values of the 
derivative of the function at the M values of the independent variable. The 
array DERI is dimensioned by variable dimension in the subroutine; therefore, 
DERI must be dimensioned in the calling program with first dimension MMAX 
and second dimension at least NCVS. 

DER2 A two-dimensional output array in which SPLDER stores the values of the 

second derivative of the function at the M values of the independent variable. 
The array DER2 is dimensioned by variable dimension in the subroutine; 
therefore, DER2 must be dimensioned in the calling program with first dimen- 
sion MMAX and second dimension at least NCVS. 

IW An input-output integer. 
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INPUT: IW is the initialization integer. Whenever a new X- or Y-array is 
input, IW must be set to -1. This condition will cause the indepen- 
dent variable array to be tested to determine if it is increasing. 
Also, certain values pertaining to the X- and Y -arrays will be com- 
puted. These values will not change unless either the X- or Y- 
arrays are replaced. 

OUTPUT: IW is an index pointer indicating that Xj^ = Xq = On the 

next call to SPLDER, the previous IW is used to begin the search 
for the interval containing the interpolation point. 

WK An array used by SPLDER as a work area. WK must be dimensioned at least 
3(NxNCVS) + 8N. This array should not be used elsewhere in the program. 

lERR An output integer error code. 

= 0 Normal return. 

= 1 The independent variable array is not increasing. A message will be 
printed, "INDEPENDENT VARIABLE ARRAY NOT INCREASING IN 
SPLDER AT POSITION UR X=XXXX,XXXX. " 

= 2 A value in the T-array is not within the limits of the X-array. There 
are no provisions for extrapolation; therefore, every Tj must satify 

Xi^TjgXn. 

Upon return to the calling program, the parameter lERR should be tested. 

Method : The method used in SPLDER is that of the reference. The reference gives 

the derivative of a matrix equation relating the second derivative of a univariate spline 
function at the given values of the independent variable to the values of the function at 
these values of the independent variable. Values of the second derivative are assumed 
to be zero at the end points. The matrix equation is tridiagonal and is solved by the 
Thomas algorithm which is equivalent to Gaussian elimination without pivoting. 
Expressions are derived for the first and second derivatives of the spline function at 
any point in an interval in terms of the values of the spline function and its second 
derivative at the end points of the interval. 

References : Greville, T. N. E., "Spline Functions, Interpolation and Numerical Quadra- 
ture," Mathematical Methods for Digital Computers, Vol. EL, pp. 156-168, John Wiley 
and Sons, 1967. 
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Subroutine PSEUDO 

Language : COMPASS 

Purpose : To create and write an appropriately named Plot Vector File. Through link- 

ages set up by an initial call to PSEUDO, all subsequent graphics data generated by the 
user will be routed through one of the PSEUDO entry points and written on the Plot 
Vector File. 

Use : CALL PSEUDO 

Subroutine INFOPLT 

Language : FORTRAN 

Purpose : To provide a one-call method of preparing plotting displays automatically. 

Use : CALL INFOPLT(IEC,N,XDATA,KX,YDATA,KY,XMIN,XMAX,YMIN,YMAX, 

PCTPTS,NXMC ,XM,NYMC ,YM,ISYM, SX,SY,XOFF,YOFF ) , where 

lEC The code for terminating the frame. 

0 (1) Used to initialize a frame and plot first curve of multiple curves per 

frame. This value of lEC leaves frame incomplete and expects addi- 
tional curves. 

OR 

(2) Used to plot second to nth curves on the frame. This value of lEC 
guarantees that all curves will be plotted in the same frame, although 
the scales for second to nth curve could be different dependii^ on 
other options selected. If the scales are different from labeled 
scales, the rescaled values of the origin and the scale increment for 
the call are printed out in the right vertical margin of the frame. 

1 (1) Used to plot one curve per frame. This value of lEC completes the 

frame. 

OR 

(2) Used to plot the last curve for the frame. This value of lEC guaran- 
tees that all curves will be plotted on the same frame, although the 
scales for the last curve could be different, depending on the other 
options selected. If the scales are different from label scales, the 
rescaled values of origin and scale increment for the call are printed 
out in the right vertical margin of the frame. 
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OR 

(3) Used to complete frame only. Generates no plotting. Special call 
with two parameters only. 

CALL INFOPLT(1,0) 

-1 Used to plot a curve with scales from previous calls on the same frame. 
This value of lEC leaves frame incomplete. If the curve will not fit on 
the plot, a new frame is set up for the remaining curves. 

Example : 

Same Plot (Scales could change) 

CALL INFOPLT(0,. . . 

CALL INFOPLT(0,. . . 

CALL INFOPLT(l,. . . 

Same Plot if Possible (If scales change, create new frame) 

CALL INFOPLT(0,. . . 

CALLINFOPLT(-l,. . . 

CALL INFOPLT(-l,. . . 

CALL INFOPLT(1,0) 

N The number of points to be plotted. 

XDATA The name of the array containing the floating-point values of X to be plotted. 

If all data within the array are to be plotted, the parameter would be expressed 
as XDATA. If it is desired to plot only a portion of the array, the desired 
beginning location is specified; for example, XDATA(4) would begin the plot- 
ting at the fourth element of the array. 

KX The interleave factor which specifies the sequence in which X-data are stored. 

= 1 Indicates that values are stored sequentially. 

= 2 Indicates that values are stored in every other location in the array, etc. 
YDATA The name of the array containing the floating-point values of Y to be plotted. 

KY The interleave factor which specifies the sequence in which Y-data are stored. 

XMIN The minimum value for X. 

XMAX The maximum value for X, 

YMIN The minimum value for Y. 

YMAX The maximum value for Y. 
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The routine checks for the first call only to determine if either 

(XMAX - XMIN) or (YMAX - YMIN) is equal to zero. When either is zero, 

the routine will scan the X- and/or Y-array to determine the limits. For 

multiple curves per display, the limits must be specified on the first call to 

include all curves since the limits from the first call will be used for all 

curves. 

If any data falls outside the limits, it will be eliminated; but a count will be 
kept of all points dropped. 

At the completion of a particular curve, if the percentage of points dropped 
exceeds the value of PCTPTS, the current curve data will be automatically 
rescaled and replotted. For multiple curves, the first curves would not be 
contained on the rescaled plot. 

If values are given for XMIN, XMAX, YMIN, and YMAX different than on the 
first call and lEC = 0, the scales could be different on the second to nth call. 
If the scales are different, the new origin and scale increment will be -printed 
out. The XMIN, XMAX, YMIN, and YMAX values are ignored on the second 
to nth call if lEC = -1. 

PCTPTS The percentage of points in any one curve that may be off scale without auto- 
matic replotting. The value of PCTPTS is a floating-point number ranging 
from 0.0 to 1. 

= 0.0 No points are allowed off scale. The data will be rescaled if any points 
are off scale. 

= 0.5 50 percent of the points are allowed off scale before the curve will be 
rescaled. Any percentage may be expressed. 

= 1.0 No rescaling will be done. 

NXMC The number of characters for the X- label including embedded blanks. The 
number of characters is calculated in two ways depending on the option 
selected in XM. 

(1) XM is written in form nHxxx. . .. 

NXMC = n 

(2) XM is the beginning storage location containit^ alphanumeric information. 

NXMC = Number of words in array multiplied by 10. Each alphanu- 
meric word contains 10 characters including blanks. 

The sign of NXMC is used to control generation of the X axis. 

>0 An axis will be drawn and annotated. 
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XM 


NYMC 

YM 

ISYM 


SX 

SY 

XOFF 


YOFF 


Language: 


Purpose: 


Use : 

X,Y 


= 0 No axis will be drawn; therefore, no annotation. 

< 0 One -inch grid lines will be drawn in addition to the axis. 

The label for horizontal annotation. This label may be expressed in two ways 

(1) The string of alphanumeric characters for the label may be written in the 
form: 

nHxxx. . .. 

(the same way an alpha message is written using FORTRAN format 
statements). 

(2) The beginning storage location of an array containing alphanumeric infor- 
mation may be used. 

The number of characters for the Y-label. (See explanation under NXMC.) 

The label for vertical annotation. (See explanation under XM.) 

An integer code specifying the symbol or mode to be used in plotting the data 
values. 

0 Draws a line (no symbols) 

The length of the X-axis in floating-point inches. The default SX is 10 inches. 

The length of Y-axis in floating-point inches. The default SY is 10 inches. 

The offset for Y-axis from frame origin. A nonzero XOFF allows room for 
the Y-axis annotation. The default XOFF is 0.75 inches. If preprinted grid 
paper is used, the user may specify XOFF = 1.0 or change offset with post- 
processor option. 

The offset for the X-axis from the frame origin. A nonzero YOFF allows 
room for the X-axis annotation. The default YOFF is 0.5 inches. If pre- 
printed grid paper is used, the user may specify YOFF = 1.0 or change 
offset with postprocessor option. 

Subroutine CALPLT 

FORTRAN 

To move the plotter pen to a new location with pen up or down. 

CALL CALPLT(X,Y,IPEN), where 

Are the floating-point values for pen movement. 
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IPEN = 2 Pen down 
= 3 Pen up 

Negative IPEN will assign X = 0, Y = 0 as the location of the 
pen after moving the X,Y (create a new reference point). 

= 999 Writes a terminating block address of 999 to terminate the Plot 
Vector File and all further processing is skipped. 


CALL CALPLT(0. 0,0.0, 999) 
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TABLE I.- ROTOR THICKNESS NOISE PROGRAM INPUT VARIABLES 


Variable 

Dimension 

Type 

Units 

BLNTN 

♦ 

Real 


CH 

♦ 

Real 

meters 

COEFFS 

♦ 

Real 


DELTET 

♦ 

Real 

degrees 

DTAUM 

♦ 

Real 


ETAMAX 

* 

Real 

meters 

JFLG 

* 

Integer 


KFLG 

* 

Integer 


MORDATA 

♦ 

Logical 

** 

N 

♦ 

Integer 


NB LADES 

* 

Integer 


NS 

♦ 

Integer 

** 

OMEGA 

* 

Real 

rev/min 

PERDM 

♦ 

Real 

** 

R 

* 

Real 

meters 

RO 

♦ 

Real 

meters 

SNDSPD 

♦ 

Real 

meters/sec 

TAUINT 

* 

Real 

seconds 

THKRAT 

♦ 

Real 


XO 

3 

Real 

meters 


Single variable, not an array. 
Not applicable, unitless. 
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Figure 4.- Definition of sweep regions. 
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Figure 5.- Rotor Thickness Noise (RTN) subprogram hierarchy. 






















tPOTOR 


xo 

= 

•25E803, 0.3, C.O, 

NBLflOES 

= 

2, 

9 

= 

.671F».Ji, 

CM 


.7E«-e0, 

OMEGA 

= 

.48E*C3, 

RO 

= 

•61E*03, 

THKRAT 

= 

0 

1 

UJ 

CO 

• 

OELTET 

= 

.4E+00, 

N 

= 

6, 

Nq 

= 

2C, 

K'^LG 

= 

r. 

f TAMAX 

= 

0.0. 

BLNTN 

= 

.1E*-01, 

JFLG 

= 

G, - 

COEFES 

= 

0.0, ,2E»01, -,2E*fll, 

TAUTNT 

= 

0.0, 

PEPDM 

= 

,106E*01, 

OTAUM 

* 

.2E801, 

SNOSPn 

s 

.34EtC3, 

MORPATA 

= 

F, 

SEND 




1 f. . 755bc3247706<.19£*00 . 36 8274 1 862 8qi26E- 0 2 

2 . 277777777777778E-C3 . 755 29 31520722C 9E+ Cli . 3764 88 34 935 7 26 2E- 02 

7 .555555555555556E-03 . 755579312092344E* Ot . 41 2 77 4 73 99568 06E-0 2 

4 . 833333333333332E-03 . 75 586 1 7262 57422E+ 0 0 . 42 330 97 361445 7 OE- 0 2 

5 . 11111111111111CE-C2 . 75614G393732565E+CC . 46 2362 00 779135 6E- 02 

6 . 138888888888888E-02 . 75641531 435 75 86E+ 0 C . 46568 70 4 7375 71 OE- 32 


7 . 166666666666665E-P2 . 756686488646952E+0C . 50 0 90 0 7260 7765 8E-0 2 

8 . 194444444444442E-C2 . 75695 391 7789762E+ 0 C . 5 041 0 1 71 45741 376- 0 2 

9 . 222222222222219E-C2 . 7572 1760 3649526E+ 0 C . 544945062 8740 86E- 02 

ir . 249999999999996E-02 . 7574775487640 0 !£♦ 00 . 53441 302 0 731 8 1 IE- 0 2 

11 . 277777777777773E-C2 . 75 77 33756344756E+ 0 0 . 565431454752224E-02 

12 . 305555555555549E-D2 . 75 798 6233 276764E+ 0 C . 559454 61 7899044E- 0 2 

13 . 333333333333326E-02 . 7582349 7511 788CE+ 0 D . 60 C 26 0 255484447E- 02 

14 . 361111111111102E-02 . 75 847 9996C982 31E» 0 C . 60 2330 16 2628 115E-0 2 

15 . 3fl8888S«8888879E-C2 . 75 872 129911945 lE+OC . 65 2646 0 1 862692 3E- 02 

16 .416666666666654E-0Z . 75 895 8890 753917E* OC .6613876642367406-02 

17 . 444444444444431E-02 . 759192778243825E* 00 . 678 877 48 31 89216E- 0 2 

18 . 472222222222207E-02 . 75942296950C225E* OC . 72662 8 27 077472 3E- 0 2 

19 .499999999999984E-0Z . 75964947310 1881E+ 0 0 . 71 289 206 8665 186E- C 2 

2C . 52777777777776CE-0Z . 7598722982941 45E+ 0 G . 76286 2 80 86 2 862 CE- 0 2 

21 . 555555555555537E-02 . 760 0 91454987670 E* 0 C . 78 54 7 30 946 77686E- 0 2 

22 .583333333333313E-02 . 760 30 6953757C 27E+ 03 . 81 3 62 928 54 8681 7E- 02 

23 .6111111111110O0E-32 . 76051 88C5839279E* 00 . 874972B74050439E-C2 

24 . 63H88888S888S67E-C2 . 76C 7 27C 23132397E* 00 , 85 5715 23 3570 625E- C 2 

25 . 666666665666643E-C2 . 760 93 16 18193645E+ 00 . 93 5 31 760 51 4966 OE- 02 

26 .694444444444420E-02 . 7611 326C4237847E+ 00 . 96 1866 129864 0 0 3E-0 2 

27 . 722222222222196E-02 . 761329995135512E* 00 . 975679 08 41 75196E- Q 2 

Figure 7.- Example of RTN output. 
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output. 
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